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IS  ABSTRACT* 


Through  various  fluctuation-dissipation  theorems,  it  is  possible 
to  relate  critical  transport  coefficients  to  time  and  space 
integrals  of  correlation  functions  of  the  relevant  fluxes.  The 
anomaly  in  the  critical  transport  properties  is  thus  attributed  to 
the  cooperation  of  two  effects:  (1)  the  anomalous  increase  in  the 
large-scale  fluctuations  of  concentration  variables  contained  in  the 
flux,  and  (2)  the  anomalous  increase  in  the  lifetimes  associated 
with  these  fluctuations  .'(-'The  fluctuation  dissipation  formula  for 
the  viscosity  in  terms  of  the  stress  tensor  fluctuations  reproduces 
the  logarithmic  temperature  divergence  of  the  shear  viscosity.  The 
excess  viscosity,  therefore,  tends  to  reduce  the  critical  Rayleigh 
linewidth  predicted  by  Kawasaki.  The  effect  of  vertex  corrections 
was  investigated  for  both  the  she.-.r  viscosity  and  Rayleigh  linewidth 
and  found  to  be  negligible  in  both  cases. 
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I.  INTRODUCTION 


A.  Phenomenology 

The  critical  temperature  of  a  binary  liquid  (a  twc- 
component  system)  is  the  temperature  above  which  the  liquid 
exists  in  equilibrium  as  a  homogeneous  mixture  of  the  two 
components  in  any  specified  concentration  .  When  the  liquid 
having  the  critical  concentration  of  each  component  is  cooled 
below  the  critical  temperature  it  begins  to  form  two  phases, 
separated  by  a  well  defined  meniscus. 

The  equilibrium  state  t  f  the  two  phases  can  be  described 
by  the  temperature,  pressure,  and  chemical  potential  of  each 
component.  However,  it  is  more  convenient  to  take  the 
pressure  P,  temperature  T,  3nd  concentration  c,  as  independent 
variables.  Near  the  critical  point  the  equilibrium  cirve 
has  the  form  shown  in  Fig.  1.  Points  within  the  shaded  region 
represent  those  states  in  which  phase  separation  takes  place; 
the  concentrations  of  the  two  phases  are  determined  by  the 
intersections  of  the  curve  with  the  horizontal  line  represent¬ 
ing  the  temperature  or  pressure  of  the  system.  Points  outside 
the  shaded  region  represent  states  in  which  the  system  exists 
as  a  single  phase  in  equilibrium.  Near  the  critical  point 
there  are  states  for  which  the  two  phases  can  be  in  equilib¬ 
rium  with  arbitrarily  close  values  of  the  concentration,  c 
and  c  +  dc.  If  y  is  the  chemical  potential  of  one  of  the 
components,  then  the  equilibrium  condition  demands  that 
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(fH T  =  °- 

r«  i  'c 


(1) 


It  can  be  shown  that  the  equal  time  autocorrelation 
function  of  the  concentration  fluctuations  is  proportional 
-  <&>  .  Therefore,  the  concentration  fluctuations  tend 

to  have  very  large  amplitudes  near  the  critical  point.  If 
we  consider  the  concentration  as  a  local  variable,  which  can 
vary  in  space  and  time,  then  it  is  convenient  to  define  the 
order  parameter  s(r,t),  which  in  this  case  is  the  concentra¬ 
tion  fluctuation  of  one  component  about  the  equilibrium 
value: 


S^t)  •£  ccr/t)  -  <C> 


(2) 


where  c(?,t)  is  the  concentration  at  space-time  point  (r’.t), 
and  <c>  is  the  equilibrium,  or  ensemble  average  value  of  the 
concentration. 

Light  scattering  measurements  on  binary  liquids  exhibit 

2 

the  phenomena  kaown  as  critical  opalescence  ,  and  critical 
3 

slowing  down  .  Critical  opalescence  refers  to  the  abnormally 
large  amount  of  light  scattering  that  occurs  in  the  critical 
region.  This  is  caused  by  the  large  regions  of  inhomogeneity 
that  exist  in  the  liquid.  The  size  of  these  regions  is 
characterized  by  the  correlation  length  £ .  At  temperatures 
far  from  the  critical  temperature  T£ ,  the  correlation  length 
is  very  small,  and  of  the  order  of  magnitude  of  the  molecular 
sizes.  However,  near  T£ ,  the  correlation  length  becomes  very 


2 


NOLTR  72-208 


large,  typically  in  the  range  of  wavelengths  of  visible 
light.  These  large  regions  of  inhoaogeneity  then  act  as 
scattering  centers  for  incoming  plane  waves  of  light,  and 
the  liquid  acquires  a  bluish  tinge,  characteristic  of 
Rayleigh  scattering,  and  then  becomes  murky,  as  it  starts  to 
separate  into  two  phases. 

Critical  slowing  down  refers  to  the  abnormally  long  life¬ 
times  of  the  concentration  fluctuations.  This  phenomenon 
is  reflected  by  the  vanishing  of  the  particle  diffusion 
coefficient  at  the  critical  point,  as  will  be  shown  later. 
Thus,  if  an  Inhomogeneity  arises,  there  will  be  almost  no 
tendency  for  the  molecules  to  move  to  the  region  of  low 
concentration,  and  so,  it  will  take  longer  for  the  system  to 
return  to  equilibrium  near  the  critical  point. 

We  can,  therefore,  think  of  the  critical  region  as  a 
region  where  the  order  parameter  fluctuations  not  only  grow 
large  in  amplitude,  but  also  become  very  long  ranged,  and 
very  long  lived.  We  will  study  how  these  phenomena  affect 
various  transport  properties  of  the  fluid. 

B.  Mean  Field  Theory 

4  5 

Mean  field  theory  *  presumes  that  the  free  energy 
density  7  can  be  expanded  in  a  power  series  in  the  order 
parameter  s(r,t).  The  first  order  term  will  be  -ys  where  y  is 
the  field  conjugate  to  s.  In  the  case  of  the  binary  liquid 
y  is  the  difference  in  the  chemical  potential  per  unit  mass 
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of  the  two  components,  that  is 


_ 

Z4  =  m,  i*)z 


(3) 


where  m^  is  the  molecular  mass  of  the  1th  component,  and 

Is  Its  chemical  potential.  In  the  absence  of  a  conjugate 
field  (i.e.,  y*o),  we  presume  the  free  energy  density  to  be 
symmetric  In  s,  and  write 


®F  =  +  ~r$*-  Z*s 

local  2 


(4) 


Minimizing  the  free  energy  density  gives  the  equilibrium 
equation  of  state 


/*-  -  as  +  bsa  , 

and  the  susceptibility  x: 


X2 


£i.m 

S-*o 


(5) 


(6) 


We  know  that  at  the  critical  point,  the  susceptibility  becomes 
Infinite;  therefore  we  presume  that  the  parameter  a  Is  temper¬ 
ature  dependent,  and  positive  for  T>Tc,  and  negative  for  T<T£ . 
We  will,  however  presume  that  the  parameter  b  has  negligible 
temperature  variation.  Therefore,  If  yi"o ,  the  quadratic  term 
In  Eq.  (4)  Is  stable  If  T>Tc,  and  becomes  unstable,  correspond¬ 
ing  to  phase  separation,  for  T<T  . , 


as  shown  In  Fig.  2. 
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For  an  inhomogeneous  liquid,  however,  we  must  also  consider 

non-local  terms  involving  gradients  of  the  concentration. 

Since  the  fluid  is  isotropic,  the  first  derivative  enters  as 

2 

a  scalar  term  proportional  to  (Vs)  .  A  term  proportional  to 
2 

V  s  transforms  to  an  insignificant  surface  term  when  integrated 

2 

over  the  volume.  Similarly,  sV  s  transforms  into  the  integral 
2 

of  (Vs)  .  Thus,  we  add  to  the  free  energy  density  the  term 

Z-1  2 

— —  (Vs)  ,  where  Z  is  a  positive  constant.  Therefore, 

*  ■  Z.*  + 


■+ 


(8) 


We  ignore  all  higher  order  terms. 

All  the  thermodynamic  properties  of  the  system  can  be 
deduced  from  the  partition  function 


where  ft  is  the  volume  of  the  system,  and 
units  in  which  the  Boltzmann  constant  is 


(9) 

8*T  ^  in  a  system  of 
set  equal  to  unity 


At  this  point,  it  is  convenient  to  expand  the  order 
parameter  in  a  Fourier  series  over  all  possible  wavenumbers, 
as  follows: 
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s<rtt )  -  Lu  ]F  S  U)  e 


irr 


Then*  ~f  y-o ,  and  ve  ignore  the  s  term. 


S*  fr'  J 


.3  a*r~t?yr 


r  r 


-*!(«♦*’«*)  IV*, 

r 


and  the  Boltzmann  factor  becomes 

-pjldr  -£l(a.,eT)IS  I2 

e  --  e  #* 


-ita.miSpI' 

e 


Thus,  by  the  equipar tition  theorem  we  see  that  the  fluctua¬ 
tions  of  different  wavenumbers  are  statistically  independent* 
and  that 


<iVl  >  =-T-rr  > 

K  K  -+-R 


where  we  have  defined  the  inverse  correlation  length 


-(aZ)  .  <14> 

The  vanishing  of  a  at  the  critical  point  thus  implies  an 
infinite  correlation  length.  If  we  define  the  spatial  correlation 


function 


)  -  <s(^>s(rD>  , 
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thea  translational  invariance  derands  that 


G&f.)  =  , 


and,  thetefore. 


G(fv.)  =  j 


h?  fce 


where 


r*.  3  r2~r* 


If  we  substitute  Eq.  (10)  into  (15)  we  obtain 


(16) 


&<?.)  =1SIL  <s  s  >e' 

r  z 


iVZ+itr'-n 


a^.<iV|2>  e 


Making  the  transformation 


a 


i 


A 

<2^ 


we  can  identify 


3r  ^ 1  sr J  ^ 


Using  Eq.  (13),  a  special  case  of  Eq.  (19)  is  just  the 
fluctuction  -  susceptibility  theorem: 

z 


\ 

~T 


1c’  T  j^,<Scr7)5(rT)>  = 


- 1 


-  a 


(17) 


(18) 


(19) 
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or 


^  ~  T  *(SVi)S(0^ 

This  suggests  the  obvious  generalization 

=  T^ft-  =  "r]A.  et<V*  <stc)stO>  . 

Finally,  if  we  substitute  Eq.  (13)  into  Eq.  (16),  then 


(20) 


(21) 


G(r)  =  (A  2L  ,  1L 

J  Cz*)3  IfV  fe1  -hit  r 


-  rcr 


(22) 


> 


which  is  just  a  screened  Coulomb  type  potential  with  a 

screening  length  equal  to  k  the  correlation  length  of  our 

syctem.  Thus  ve  see  that  in  the  critical  region  where  the 

correlacion  length  becomes  very  large,  the  range  of  the 

correlation  function  tends  to  infinity.  Eq .  (22)  is  a  result 

first  arrived  at  by  Ornstein  and  Zernicke^  in  1914, 

He  recall  that  the  above  results  are  only  valid  when  we 

can  neglect  the  fouxth  order  term  in  the  free  energy  density. 

However,  as  we  get  closer  to  the  critical  temperature,  the 

coefficient  of  the  quadratic  term  will  tend  to  zero  and  the 

fourth  order  term  will  start  to  dominate.  He  observe  that  a 

change  of  scale  by  a  factor  X  in  the  order  parameter  Induces 

4  -1  -1  2 

the  transformations  b-*bX  and  Z  +Z  X  ,  so  that  the  quantity 
2 

bZ  ,  which  has  units  of  reciprocal  energy  times  length, 
remains  invariant.  Thus  if  ve  define  a  length 

t 


K  -  , 


(23) 


then  we  expect  that  the  above  results  of  mean  field  theory 
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should  hold  when  t^-<<1, 
larger  than  £0* 


and  start  to  break  down  when  £  becomes 


C.  Hydrodynamics 

It  Is  necessary  to  employ  hydrodynamic  equations  of  mo¬ 
tion  in  calculating  the  time  dependent  correlation  functions^ 

of  concentration  and  velocity.  TLe  time  dependence  of  the 
3 

concentration  is  derived  from  the  continuity  equation 

—  \f-f  -  O  (24) 

-dt  J  ' 

where  J  is  the  particle  current,  given  by 

J(r*,t)  -  Slr’tJvC^t)  ,  (25) 

where  v(iT,  t)  is  the  fluid  velocity  at  space  time  point  (r,t). 

It  is  also  possible  to  write  a  constitutive  equation,  analogous 
to  Ohm’s  Law  in  electricity: 


f  -  -XV/.  ,  <“> 

where  X  is  the  particle  conductivity.  If  we  make  use  of  the 
definition  of  the  susceptibility,  then  we  can  write 

-  ^T)  =  \  ^  .  (27) 

and  Eq.  (24)  becomes 

—  -  j)  (28) 

* 

where 

D  «  (29) 
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This  is  just  the  diffusion  equation,  with  D  equal  to  the 
particle  diffusion  coefficient.  More  generally,  we  have  the 
equation  for  the  wavenumber  dependent  concentration  fluctuation 


it  sr‘tJ  -  - 


(30) 


where  r^» is  the  decay  rate  of  the  kth  mode,  given  by 


2 

x  k 


(31) 


The  solution  of  Eq.  (30)  is  just 


-  r  t 

s^a)  -  e  s£,(o)  , 


(32) 


and 


3r 


It,.)  "  <s..^)5t,(s)>  =  <isj'>e 


(33) 


To  determine  the  velocity  autocorrelation  function,  one 
makes  use  of  the  equipar tition  theorem  to  calculate  the  equal 
time  value,  and  the  Navler-Stokes  equation  to  calculate  the 
time  dependence.  First,  we  Fourier  analyze  the  velocity  field, 

itr 


(34) 


and  resolve  v^-(t)  into  orthogonal  components 


v_<t)  -  I  f  y ■  n) 

K  JU-t 


(35) 


We  have  only  two  directions  of  polarization,  since  we  will  only 


be  interested  in  the  transverse  velocity  modes,  with  e 


ky 


representing  a  unit  vector  along  the  direction  of  polarization, 
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such  that 


*V°  - 


u  -  i,  1  • 


Ve  can  vrice  the  kinetic  energy  of  the  fluid  as 


20*7,  ]p.<v«-  V*-  V  W  J1 


ar-r^r* 


»ifizivr/. 

!  a-/  K 


The  equlpartition  theorem  then  implies  that 


J<'Vr/>  =  T  - 


<'VJI>'=  T 


where  p  is  the  mass  density  of  the  fluid.  Again  we  have  the 
situation  that  different  velocity  modes  are  statistically 
independent.  To  determine  the  time  dependence  of  the  velocity 
modes  we  make  use  of  the  linearized  Navier-Stokes  equation: 


v  rovo. 


where  T  is  the  stress  tensor,  whose  cartesian  components  are 
given  by 


11 
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where  p  la  the  hydrostatic  pressure,  and  n  and  C  are  the 
shear  and  bulk  viscosities,  respectively.  If  we  separate  out 
the  transverse  components  of  the  velocity,  that  is  the  part 


for  which 


V  -  7  c^t)  ^  o  s 


then  Eq.  (39)  becomes 


—  v.^t) 

Bt  4 


5  V°  = 


lfe2v  to 
P 


and  thus  we  can  write  the  velocity  correlation  function  as 

T  &2  jt  l 

<Vt«’Vt,)>“  p e  '  •  <44: 

He  note  that  the  velocity  equal  time  correlation  function  is 
not  k-dependent,  whereas  the  equal  time  concentration  correla¬ 
tion  function  has  the  k-dependent,  and  algo  temperature  depend¬ 
ent  0rn8tein-Zernicke  form.  We  also  note  that  the  velocity 
modes  will  decay  much  faster  than  the  concentration  modes  in 
the  critical  region,  because  we  expect  that  D  will  go  to  zero 
strongly  while  ^  will  either  have  a  weak  divergence  or  stay 
finite. 

In  this  chapter  we  have  developed  the  tools  of  mean  field 

12 
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theory  and  hydrodynaaics ,  which  will  be  indispensable  in  our 

analysis  of  the  critical  variation  of  transport  coefficients. 

In  each  case  we  will  relate  the  transport  coefficient  to  an 

integral  over  all  space  and  time  of  a  correlation  function  of 

8  9 

the  relevant  fluctuating  quantities  '  .  To  evaluate  these 
integrals  we  will  need  to  use  the  concentration  and  velocity 
autocorrelation  functions,  which  we  have  just  derived. 
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II.  CRITICAL  SHEAR  VISCOSITY 


A.  Fluctuation-Dissipation  Formula 

The  wavenumber  and  frequency  dependent  shear  viscosity 

can  be  computed  from  a  fluctuation-dissipation  formula 

relating  the  viscosity  to  the  Fourier  transform  over  space  and 

time  of  the  autocorrelation  function  of  the  local,  time 

dependent,  off-diagonal  component  of  the  stress  tensor 

3”  (r*, t)10,1*.  That  is 

*  2  ’ 


If  *  -£=  j d\,  e 


(1) 


where  T  is  the  absolute  temperature  in  a  system  of  units 

in  which  the  Boltzmann  constant  is  set  equal  to  unity  (k_El) , 

1) 

and  where  r^-r^-r^  and  t21*t2”tl  are  relative  space  and  time 
separations. 

To  calculate  the  excess  viscosity  in  the  critical  region 
it  is  necessary  to  relate  the  local  time  dependent  stress  tensor 
to  the  order  parameter  fluctuations,  which  become  very  large 
in  amplitude,  extremely  far  ranging,  and  very  long  lived. 

This  is  accomplished  by  identifying  the  Ginzburg-Landau  free 
energy  density  with  a  Lagrangian  density,  and  using  a 
canonical  expression  relating  the  stress  tensor  to  the 
Lagrangian  density. 

The  most  probable  configuration  of  a  system  is  that 
which  maximizes  the  partition  function  Z^,  or  equivalently. 
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minimizes  the  total  free  energy 


3 

r  , 


since 


(2) 


where  £  denotes  the  sum  over  all  possible  configurations  in 
,Sj 

space  of  the  order  parameter  s.  Thus,  if  we  consider  the  free 
energy  density  $  to  be  a  functional  of  the  order  parameter 
and  its  gradient,  i.e. 


T  -  Hs,s  •] 


i  -  ij  2,  ^  , 


(3) 


where 


(4) 


then  the  problem  of  the  most  probable  configuration  of  the 
system  is  equivalent  to  minimizing  the  free  energy,  and,  thereby, 
becomes  a  problem  in  the  calculus  of  variations,  where  the  vari¬ 
ation  of  the  integral  over  all  space  of  ?  is  zero: 


fji’r  TU.SJ  —  O  , 


(5) 


or,  the  functional  derivative  of  T  with  respect  to  s  vanishes: 

s? 


cTs 


.=  o 


Eq.  (5)  then  becomes 

S  [dV  J<*V f65;]  -  j  j 
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If  we  cote  that 


*\  -  -  £<**>  , 

and  if  we  integrate  the  second  term  in  Eq .  (7)  by  parts. 


ignoring  surface  terms,  we  get 


=  o 


For  arbitrary  variations  6s,  this  leads  to  the  Euler- 
Lagrange  equation 


3?  <r-  c>  7)5  r. 

«  -  -  °  •  (10) 

Therefore,  /e  can  employ  a  canonical  formalism,  with  the 

free  energy  density  analogous  to  the  Lagrangian  density  in 

12 

scalar  field  theory  . 

We  can,  therefore,  use  the  absence  of  explicit  coordinate 
dependence  of  the  free  energy  density  IF  to  derive  an  expression 
for  the  stress  tensor  .  Under  an  infinitesimal  displacement 


X;  =  X-  6t 

the  free  energy  density  changes  by  the  amount 

y*  -  r-r  -  Z  t;  57..  - 


If.  however,  IF  has  no  explicit  coordinate  dependence,  then 


?£  r.,  .  y  S  i\. 


w-  -  Z 


2>S:  •  J 
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where 


3  .J 

Ss  =  S<T riT)  -  Mr  )  "  ax-  ~  ^4  Si  • 

j=i  1  f' 


But,  by  the  Euler-Lagr ange  equation,  Eq .  (13)  becomes 

$$  “  £ Lfo  if,  )^<j  ■**  (?x4-  *  S  )J 

j*« <r' 

If  we  equate  Eqs.  (12)  and  (15)  we  obtain 

1 t '  A 1 S.  sr  V r  1  •  5  • 

or 

I  &.  T  •  =  O  ,  <17 

i-l  J 

where  we  identify  with  the  stress  tensor,  which  we  define 


T  ^  ^  T 

J;a  ”  ast  j  *j  • 


In  particular,  let  us  now  consider  the  Ginzburg-Landau 
equation  for  the  free  energy  density 

=  fs!4-  %-<**>*-  ~5--/‘s  j 
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where,  in  the  binary  liquid  we  identify  the  order  parameter 
a  with  the  excess  local  concentration,  and  y  with  the  chemical 
potential,  then  the  Euler-Lagrange  equation  provides  us  with 
the  "equation  of  motion,"  or  in  this  case  the  equation  of 
state : 

u  -  as  *-  faS3  -  Z  ,  (20) 

and  the  off  diagonal  component  of  the  stress  tensor: 


o;2cr,t)  = 


7>T  ^ 


(21) 


We  note  the  ys  term  in  the  diagonal  components  of  the  stress 
tensor,  so  that  the  divergence  of  the  stress  tensor  will  contain 
a  term  V(sy),  which  is  recognizable  as  the  force  per  unit 
volume  caused  by  the  gradient  in  the  chemical  potential  density. 
This,  in  turn  is  equal  to  the  time  rate  of  change  of  momentum 
density,  or  momentum  density  flow,  in  a  region  of  fluid  Thus, 
we  can  justify  the  choice  of  unity  as  the  arbitrary  constant 
factor  involved  in  the  definition  of  the  stress  tensor,  Eqs.  (16)- 
(18). 

From  now  on  it  will  be  convenient  to  work  in  wavenumber 
space,  rather  than  configuration  space.  This  is  accomplished 
by  Fourier  transforming  the  order  parameter  as  follows: 
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—  rv/z  y  ^(tj  g 


;r-F" 


Vt;  n 


r  3  ^  -t£r 

Vz  J dr  S^t)e 


where  SI  Is  the  volume  of  the  system.  In  the  thermodynamic 
limit  It  Is  possible  to  transform  the  sum  over  IT  Into  an 
Integral  over  lT  by  the  prescription: 


<rl 


J  (**)* 


where  the  upper  limit  of  Integration  is  a  Debye  cutoff. 

Similarly,  if  we  Fourier  transform  the  local  stress  tensor: 


■^2  if  6 

I? 


>k-r 


rj(i)  ■=  jpkJaV  J^rt)e 


-l&r 


then  by  Eqs.  (21)  and  (22),  (26)  becomes 


r/tt)  -  jp  Z  (AkrCt)  Sect)  se_eCi).  (27) 

We  can  now  rewrite  the  fluctuation-dissipation  formula 


for  the  shear  viscosity,  Eq.  (1): 


«ot2(  | 
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Translational  invariance  of  the  stress  tensor  demands  that 
k* ■  -k ,  so  that  Eq.  (28)  becomes 

V  >  -  sK .***  ' r" 

*  It 

=  ^  ^  SKf 

.  «*»> 

To  calculate  the  anomalous  part  of  the  sh-iar  viscosity 

in  the  critical  region,  we  substitute  the  contribution  to 

JT  q(i)  of  the  concentration  fluctuations,  namely  Eq.  (27), 

X  z 

so  that  Eq.  (29)  becomes 

~  ~ zT&a  1a, e<  4‘ 

r  c 

X  <SrCt,)Sf  rCtJ)^,Ct,)S-4.(e,)>  .  (30) 

At  this  point  »e  have  to  consider  the  correlation 
function  in  the  above  equation.  We  must  correlate  two  com- 
pouents  of  the  order  parameter  at  time  t^  with  two  other 
components  at  time  t^«  We  will  make  a  decoupling  approximation, 
which  we  will  justify  later.  This  approximation  suggests  that 
we  treat  this  correlation  function  as  the  sum  of  all  possible 
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products  of  two  correlation  functions  of  the  fora 

<s.  (t,)c.  (t.)>.  In  Feynman  diagram  notation,  this  is 
*1  1  *2  1 

equivalent  to  considering  only  graphs  of  the  form: 


(31) 


Later  we  will  show  that  higher  order  graphs  do  not  contribute 
to  the  shear  viscosity.  Thus  we  will  write 


<S-.U2)  Sf  -CtJ  sr,( t.)  5 


VV  f(t2>  ?f )>  +<^.(ti)S^.)><5rr(ti)  S-tt.  )> 


+  sr-r.?  (t'^' ^ 


“  <s,-  ft.)  V°>  <sf.r  (*,)  sf^(t,  »  (J, 


If  we  assume  that  the  order  parameter  fluctuations  decay 
exponentially  with  a  wavenumber  dependent  lifetime  we  can 

write  each  time  dependent  correlation  function  as  the  equal 
time  correlation  function  times  a  decaying  exponential  as 
follows : 
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<\.i^ sr«,)>  =■  e 


f  t,( 


(33) 


If  we  substitute  Eqs.  (32)  and  (33)  into  (30)  ve  obtain 


=  -2-f^  1 Z  A.  (■!  (fc-4  Mfe-  4  )<i\.  i’>  <i  sf_  ,il> 


(34) 


In  the  above  equation  we  will  take  q  to  be  along  the  z-axis , 
so  that  Eq.  (34)  becoaes 

A*||*(«)  =  2T&CI  ^  %i  * i  )  “  ^ 


(35) 


To  exhibit  and  aake  use  of  the  synaetry  in  the  above  expression 
it  is  convenient  to  make  the  substitution 

F 


£•%-*■ 

C  -  -A* 


£ 

L. 


(36) 


then 


Ar|*(u) 

wIuiplpr  . 


X 


(37) 
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If  ve  note  that  ve  can  write 


4 ’it**'-*.)  4) 


(38) 


then  only  the  first  term  will  contribute  to  Eq.  (37),  because 
the  second  tera  will  antisymaetrize  (37).  Then 

=  zrlfa  ‘ fx( 1 V 1  ><' V 1  > <J^uy7rJj‘  .<W 

Now  we  are  in  a  position  to  apply  the  above  formula  to 
a  few  special  cases. 


B.  Temperature  Dependence  at  Zero  Frequency  and  Wavenumber 

“’j  -* 

At  2ero  frequency  and  wavenumber,  1  »  -1,  (ii»o,  and 
Eq.  (39)  becomes 

•2. 

.i1 

i 


Av(o)'  j 


(40) 


where  K  ^  is  the  correlation  length.  In  the  critical  region 
we  will  use  the  Orns tein-Zernicke  formula  for  the  wavenumber 
dependent  order  parameter  correlation  function,  namely 


1r  "  ** 


i.  T 


(41) 


**  -8-  r 


For  the  linewidth  of  the  order  parameter  fluctuations  we  will 

13  14 

use  a  modification  of  the  Kawasaki  '  formula: 


a  (d-\ t  K  i 


(42) 


where  n  is  the  non-anomalous ,  background  shear  viscosity, 
o 

which  is  the  component  of  the  shear  viscosity  not  influenced 
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by  the  critical  fluctuations.  The  Kawasaki  scaling  function 


K(z)  is  a  slowly  varying  function  with  a  minimum  value,  K(o) , 


3ir 


of  1,  and  a  maximum  value  K(»)  of  (-1.18),  explicitly  given 


by: 


KOO  =  f +  (j - •?)**»'*]  . 


(43) 


Experimentally,  it  has  been  found  that  the  linewidth  behaves 


as  if  the  scaling  function  were  constant  and  equal  to  its  high 

3tr 


wavenumber  limit,  and  n  replaced  by  an  experimentally 

o  o 

fitted  parameter  rj,  about  29%  larger  than  the  background 


viscosity  no.  Therefore,  we  assume  that 


r  =  I£ 

■t  IG^ 


Substituting  Eqs.  (41)  and  (44)  into  (40)  we  get 


_  IG  y 

n  ~  O-  4^ 


a  -  i 

We  first  calculate  the  angular  average: 

2Tt 


{ ^r- )  /=  j suede j a*  [s^e  <*s a  &s2e] 

x  9,p  *  o 


ITT 


~  jtos*©  sin*©  40  j  Co 5*4> 


so  that  Eq.  (45),  upon  transforming  the  sum  to  an  integral 
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(45) 


=  =  7(3  "?)  =  Ts  >  <46) 
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thcougb  Eq.  (24) ,  becomes 


\  JL_  f 

~ir  =  H*- W)  (K%F)s/i 


-  JL  {  l  dC 

'  isn') 


(47) 


To  evaluate  the  above  integral  we  make  the  trigonometric 
substitution: 


X  =  K  ttut  9 


(48) 


so  that 


r  iUi  =  i 

Jc «Ktyu  J 

o  © 


lari  K 


-bSe 

Sec  © 


T  - 
XidH 


-I  ■] 


ton'^'  ^  W  ■“ 

■S  f  ciG  -  [  Siu1©  C 

CoS  © 


COS  0  cl  © 


-  [-Sin  0  +  Cn  tan  v  ~  ~  ^  Siq3 &  j 


tan 


(49) 


Iu  the  critical  region,  where  the  correlation  length  becomes 
kD 

very  large,  —  >>1,  and 


« 


tan  — 


jr 

Z 


K 


(50) 
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Therefore: 


♦a  M 

(  JjL L  ~ 

l  ~ 


Jin  ~ 


to  the  lowest  order  in 


zero  wavelength,  we  have: 


Thus  for  zero  wavenumber  and 


<*?  -  JL  \  e  2ka  r.<-o  l 

n  "  i5irll>  IT  S  J  > 


with 


C,0>  ~  -%r  ~  -1*333  . 

\  d 

Assuming  that  ic  vanishes  like 


K  « 


T-Tt, 


then 


_  a  t  .  ,  ,  .  l 

--  — ;  1  -  V  JLv\\  —zr~  -+-  Constant  I 

X  >s^L  1  T<  ‘  -J  * 

Anfic)  T-Tc 

Therefore,  if  — L~  is  plotted  against  In  — - —  ,  the  slopi 

n  c 

of  the  line  that  results  should  be  — 

151T 

The  critical  shear  viscosity  of  the  binary  system 
3-methylpentane-nltroethane  has  been  investigated  by  Stein, 
Allen,  and  Allegra^,  who  found  a  slope  of  -.039.  The 
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critical  exponent  v  for  the  divergence  of  the  correlation 
length  has  been  found  to  be  equal  to  .616  for  the  3-methylpen- 
tane-nitroe thane  system,  which  would  yield  a  theoretical  slope 
of  -.033,  a  value  15Z  lower  than  experiment. 


C.  Four ier-Laplace  Transform  of  Critical  Shear  Viscosity 

In  the  calculation  of  the  correction  to  the  Rayleigh  line- 
width  due  to  the  anomalous  shear  viscosity  it  will  be  useful 
to  know  the  Fourier-Laplace  transform  of  the  critical  shear 
viscosity,  defined  by 


o 


This  will  only  modify  the  frequency  factor  of  Eq.  (39)  as 
follows : 


-  j 


<isri2><!  V1') 

r  rv  *  rr 


(57) 


We  will  calculate  this  quantity  at  the  critical  temperature 
T-Tc,  where 


(58) 


(59) 


If  we  convert  the  Laplace  frequency  to  wavenumber  units, 
by  the  equation 


X 


X 
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Q.*(Vf 


(60) 


and  Insert  Eqs.  (24),  (58),  (59)  and  (60)  into  (57),  we  obtain 

fL 


aaV(<?)  =  (*£c  */(/;- of 

Let  us  now  consider  the  important  special  case  q*cT. 


(61) 


we  can  say 

A'— A, 

and  Eq.  (61) 

become! 

—  | 

A 

[  a3* 

i/c, 

z 

l 

a 

IT3 

t 

1  v 

i  e\  i 

Q* 

_  Z_ 

* 

4tt  ( 

**  jt 

it 

•5  J 

cK 

_  _s_  .  a- *1 

-  *w‘  3  H  (%) 


(4? 

3 


g  D 

rz  ~zzi  Jen  ~pr  (62) 

I5«4 

for  small  enough  Laplace  frequency  to  consider  Q<<kQ.  Again 

g 

we  get  the  ——7  times  a  logarithm  term,  with  the  inverse 
15TT 

correlation  length  replaced  by  the  frequency  in  wavenumber 

18  19 

units,  as  we  would  expect  from  the  principle  of  dynamic  scaling  * 
Because  of  the  convenient  definition  of  Q,  we  have  no  additive 
constant  to  the  logarithm,  as  we  had  in  the  case  of  zero  wave- 
number  and  frequency.  Equivalently  we  can  say  that  (In  k-C^0^) 

-c(°)  11 

is  to  be  replaced  by  In  Q,  or  that  Q  scales  as  Ke  i.e. 
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<«»  J 


C'"J  -  ^3 

Qe1  -  Qe  - 


Returning  now  to  the  aore 
or  graphically: 


general  case,  where 


it  is  convenient  to  transform  the  integral  over  y  ■>  cos6 , 
an  integral  over  £’  through  the  law  of  cosines: 


iMMMaNMnittiiitf 


P*  ..-i  -a 


(63) 


f 

» 


to 


(64) 


(65) 


(66) 


(67) 


(68) 


(69) 


"  ■ . . . TUI  nil—  M 
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Noting  that  the  integration  over  the  azimuthal  angle  gives  a 

r2ir  2 

factor  of  J  cos  $d$  -  ir,  and  substituting  Eqs.  (64)-(69) 


into  (61) ,  we  find : 


A  ♦  4a  (ft+W 

a  ^  o  1  ‘H^ai+£i5-+  £$*)  .  (70) 


p 


If  we  transform  to  dimensionless  units,  where 


, ..  £ 
u~  % 


p=  WqP  =  ±P 

T  7lV»«i 


then  Eq.  (70)  becomes 


4j»  (*♦») 


-  _!_  (  du  f  4a!,  g(.'%  u'O  - 1-  c«a- 

^  j  a  J  u‘  cuVu'V-p)L  J  ; 


0 


where  the  region  of  integration  is  shaded  in  Fig.  6. 

The  double  Integral  in  Eq.  (73)  is  more  readily  evaluated  by 
a  45°  rotation  of  the  axes;  letting 


X  -  dtV  u  1 
*'  *  u'-  u  ‘ 


u'- 

U’  t(x-x'jj 
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then,  the  Jacobian  of  the  transf oruation  is 

7>C*jX‘)  I 

and 

HU.'  —  IfC**-  y  (77) 

tl*2- Ul  -  (u''U)(U.'+  U)  =  XX'  ^  (78) 

j£(uW*)- i- Cu^-u^J  =  x2+x'z-  i  -  xix'z 

=  (x-«>0-x'z)  <79> 

uK  u'*  =  -g-(x3-t-  3xx'2  )  .  <80> 


Substituting  Eqs.  (75)-(80)  into  (73)  we  obtain 
*  , 

Ar)-y(p)  =  |  *  UMx  xix'i(xl-O0- x,2J 

a  J  «x*-x-o  p) 


1  -t 

24 


—  —  f ^  l  i  \[  dfx'x'* 


0  -  x'2) 


Cx*-x'2X*3+:sxx'MO 


(81) 


where  the  function  f(q,P)  is  defined  by: 


IS  /  %  , 

Z  J  4x x  1 LX"’)Jc*U-t)Cx% 3xx''t-P} 

.  O 


(82) 
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Noting  that  the  double  integral  in  the  above  equation 
ranges  over  values  of  x  larger  than  unity,  and  values  of  x' 
smaller  than  unity,  we  can  approximate  it  by  neglecting  x' 
and  unity  when  added  or  subtracted  from  x,  so  that  a  first 
approximation  to  f(q,P)  is: 


-  if) 


P 


I  +  P 


?\‘/3 


c 


(83) 


Again  we  recover  the  logarithmic  dependence.  If  either 

the  wavenumber  or  frequency  is  very  large,  such  that  q  and  Q 

are  of  different  orders  of  magnitude,  then  the  ratio  of  the 

excess  shear  viscosity  to  the  total  shear  viscosity  is 
2 

equal  to  8/15ir  times  the  natural  logarithm  of  the  ratio  of 
2kD  to  the  larger  of  q  or  Q. 

We  note  that  f(q,P)  depends  on  q  only  through  the  upper 
limit  of  the  x  integral  in  Eq.  (82).  Thus,  if  we  subtract 


f  (q , P )  from  f(q,P),  then  at  least  to  order 


2  k. 


we  can  extend 


the  limit  to  infinity,  having  subtracted  -"it  the  divergent 
part.  If  we  define 

s P)-  =  Vp>  - 
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then  the  scaling  function  0^  will  only  be  a  function  of  P, 
which  is  proportional  to  the  ratio  of  the  Laplace  frequency 
to  the  cube  of  the  wavenumber.  Therefore,  if  we  subtract 
Eq.  (83)  from  (82)  and  extend  the  integral  to  infinity,  then 
we  have 


1S.,~  ,f'.  ,  .S  I 

-  yUxx1  Ux  x'  (i-/.  ^  “  — 

Z1  i  p)  ApJ 


2  *  i4*  X  °"X  ^[(>-T(  ^(x3,3AX* 


3xx 


.1 


1-P)  (/-PXx-axy'rP 


(85) 


It  is  possible  to  perform  the  x'  integration  if  we  let 

<Ux>  ^CX3i-  P)'1  ,  (86) 

b(x)  S  (3*)'1  ,  (87) 

then 


,5-f  [  •  J  Q-*‘> 

«  A 


q2Ca2+b2x'r) 


bY*0  +  £(aVbV)xZ1 

+ bVO  +  b*  tf'2)  J 


I 


f'axVVVX'-a.xV1)  l 

J  ~<$  -  1 
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where. 


and 


(89) 


(90) 


To  evaluate  J^(x) ,  it  is  convenient  to  sake  a  hyb^rbolic 


substitution: 


X  -  X  (ink  B 


<ix  -  xseclfte  tiO  , 

so  that 

loo 

IT 

~*J  tWfeO-Zx^k*©  +X‘1U^6)d9  ,  (92) 

o 


Noting  that: 


d€>  =  JO- Scch*6 )JG>  ~  O-tiinkO  ) 

J t/}r>h* 0^0  ~  J^Amk  &(j-$ech  ®  £>  ‘~3’tftrtK©  ^ 

J  tkrih  &  ~  Jt&rtA,  ®0~  Sf’Ci i  b)<H&  -  3^«li  Q  ~  5  t(i«k  B  ^ 
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I 

we  compute 

J,  (x  >  =  6-  X*)Z(y  tiUilt ‘*5?  ”  1 )  +  *k  Or  X*  )  “ 

*  0  -  «*)*(*  tank  *X  %G~**)  +  —  ,  (94) 

j 

1 

To  evaluate  ^(x),  it  is  convenient  to  make  a  trigonometric 
substitution: 

Y  fcv»  B 

x'2  *  ^(l  +  W€>)  -  ^  Se/p  (95) 

dx'  -  %  S6C2<^d^ 


so  that 


T  /  X  -*7  "  *»?€>$*"•»(■ - xWfpJP 

JjUO  J.  ^Tsecze 


(96) 


ta«  “ 

=t jlWa  -  (.£  -  Sfa'e  - 1? x*  ta"l'e’  ] • 


Noting  that 


JiknOde  =  jcsec?€>- i)de  =•  6t»B  -  61 


Jtkn,&<l&  ~  JtA«2£>  (secG>- Odf?  -  +  \t(Ui  © 

(SfvV-O  r  C ©)-  +  s^w  £ 


(97) 
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we  compute 


T,t.>  -  W  ^ 

-  o <98> 


Therefore: 


T(x).  J,w  =  0-x,)f('-x‘X*M-Ji-')>0+  )] 

i5  V  d*  J  , 


and  by  Eq.  (88) 


*» 

•qi?)  =  -  f  |^:{o-xo{(,-x*)Olr^i.  i  -  ■) 

i 

4f)(«- (100) 


To  study  the  behavior  of  a^(P)  for  very  large  values 

of  P  It  is  simplest  to  refer  back  to  Eq.  (85),  and  note  that 

2 

for  large  enough  values  of  P,  we  can  neglect  the  term  3xx' 

3 

compared  to  x  +P,  so  that  Eq.  (85)  becomes,  approximately 


>  0 
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,-2/3 


He  can  see  that  the  above  expression  will  vary  as  P 
lowest  order  in  F-1.  In  fact  if  we  define 


u  p>-A£ 


then 


iJl 

J  6cV 


x—  ^  -  —  TCP) 
py  <JP  JA  ;  * 


Evaluating  I(P)  we  find: 


up)  *  J*b  ’  i  pfefe  * B  (* ' w  #s^ 1 


~  22-  p  ih 

3f3  r 


and 


.  2-  2  TT 

jpXtP)  -  y  57%  p  - 


A 


Thus  Eq.  (101)  becomes 


<yp) 


I 

- j  P!A  If^xVW1*^2)-  ^J^VV*') 


►-!ip2/3  2£  \ 

2  r  3^3  7  ^-7J 


Ztr  p'^3 

3/3 


-2/3 

-  1.20*1  P 


(101) 


(102) 


(103) 


(104) 


(105) 


] 


(106) 
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The  scaling  function  was  calculated  for  various  values  of 
P,  by  numerically  integrating  Eq.  (100).  In  the  process  we 
transformed  the  infinite  range  of  integration  to  a  finite  one 
by  the  transformation 


X  -  tan  4> 


(107) 


and 


«  % 

- J jowo  sec  4'  Jcfr 


(108) 


% 

The  results  are  plotted  in  Figs.  (3) -(5). 

He  note  that  for  the  special  case  P*0,  we  have 

<r(°)  7  )J 


-0.633  . 


(109) 


Eq.  (85)  was  also  used  to  calculate  the  initial  slope  of  the 
scaling  function,  as  follows: 

«*»  t  t  -j 

°it(o:)=:  -  Cxi-X'0(xiV3Ox2] 

jax  x'Vx'1)  -  T ]<** *' ^'"x^|(Tx 


- 

2 


»  o 


1<rf'  r°°  ix(x-t ) 

-  3  2.  J^x*  ^  J  <A-x,aX*V3x'iJ 

o  • 

i  -  -  *•' )  J*  1^7*  -  .Ti7‘] ' 
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-  0.25H 

Therefore,  for  small  P 


(110) 


(TIP)  ^c- 0*4,33  +  O.ZSMP  (lu> 

t 

The  integrals  in  Eqs.  (109)  and  (110)  were  evaluated 
numerically.  Details  of  their  evaluation,  along  with  the 
general  evaluation  of  o^(P),  Eq.  (100),  can  be  found  in  the 
Appendix. 

The  wavenumber  dependence  of  the  critical  shear  viscosity 
in  the  low  frequency  limit  can  now  be  written  as 


^  =  4Jfnf  -0.«3| 

^  'STT  l  %  ]  ) 

or  equivalently  we  can  say  that  K  and  q  scale  as 


~in  k  -  1.333  Cr 


0.4,33 


(112) 


or 

-0.700 

k 


K 


(113) 


In  this  chapter  we  have  calculated  the  critical  shear 
viscosity  by  use  of  the  fluctuation  dissipation  theorem,  and 
have  verified  the  experimentally  observed  logarithmic  divergence. 
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We  have  alao  shown  how  the  correlation  length,  wavenumber, 
and  frequency  scale  with  respect  to  each  other. 
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III.  WAVENUMBER  DEPENDENCE  OF  THE  CRITICAL  RAYLEIGH  L1NBWIDTH 


A.  General  Formalism 

The  Rayleigh  linewidth  is  the  decay  rate  of  the  q* 
wavenumber  mode  of  the  concentration  fluctuations.  It  is 
given  by  the  formula 


where  Dht is  the  q-dependent  diffusion  coefficient.  The 
diffusion  coefficient  is  equal  to  the  ratio  of  the  conductivity 
to  the  susceptibility  x^*  The  susceptibility  is  proportional 
to  the  equal  time  value  of  the  correlation  function  of  the  q- 
mode  of  the  order  parameter. 


*f  =  T<|Sf|1> 


(2) 


By  the  use  of  a  fluctuation-dissipation  formula,  we  can  relate 

the  conductivity  to  a  space-time  integral  of  the  correlation 

14 


function  of  the  solute  particle  current 


That  is, 


Xf  -  Zfjdt,,  j  A.  J-i' ^ ^  •  ( 


3) 


where  J^(r,t)  is  the  component  along  the  q-direction  of  the 
local  solute  particle  current,  which  is  equal  to  the  second 
order  product  of  the  concentration  fluctuation  a  times  the  local 

-to 

fluid  velocity  v, 

A 


(4) 
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Therefore,  any  simultaneous  fluctuations  In  the  concentration 
and  velocity  fields  lead  to  nonvanishing  contributions  to 
the  current,  according  to  Eq.  (4).  Such  fluctuations  In 
the  current  In  turn  give  a  contribution  to  the  conductivity  by 


Eq.  (3). 

At  this  point,  It  Is  useful  to  make  a  decoupling 
approximation,  separating  the  concentration  fluctuations  from 
the  velocity  fluctuations.  We  can  do  this  If  we  note  that  the 
transverse  velocity  fluctuations  do  not  lead  to  any  concentration 
changes  In  the  liquid,  and  as  a  result  have  much  lower  frequency 
components  than  the  longitudinal  modes.  Consequently,  the 
transverse  velocity  modes  make  a  much  greater  contribution  to 
the  Integral  In  Eq.  (3).  Therefore,  we  can  Ignore  the  longi¬ 
tudinal  velocity  modes,  and  decouple  the  transverse  inodes  from 
the  concentration  modes.  This  assumption  now  lets  us  write 
the  Integral  In  Eq.  (3)  In  factorized  form,  as  follows: 

cr> 


It  will  be  convenient  for  us  to  work  in  wavenumber  and 
frequency  space.  Therefore,  we  will  Fourier  transform  the 
concentration  and  velocity  fields  as  follows: 


Sir,  t )  Q'/i  2.  |  i~  e 


iCF’f*- v t) 


\  W  V")e 
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If  ve  sake  use  of  the  translational  Invariance  in  space  and 
tine  of  both  the  concentration  and  velocity  inodes: 

<  S ^»)  5_0O>  =  2nS(»+ v‘)i  <IS 

a  t  *  c  * 


and  we  recall  the  Fourier  transforms  of  the  delta  functions 


=  25"  8  (to*  v)  , 


ur^-p-r 


=  QJl 


T*  * 


r.pe 


then  Eq.  (5)  becomes 


*r 

For  the  concentration  fluctuations  we  will  use  the 
Ornstein-Zernicke  form,  with  a  decay  time  of  f.  That  is, 


<iy^)|1>  =  J  <rt2i  e‘“’t'{sr(t‘)Vt'  J> 


■0V,>&^e!i’W 
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where  <Js£,  \>  is  the  equal  time  correlation  function,  and 
j£(w»r;?i)  is  the  Lorentzian  function  of  the  frequency  w  with 


a  half-width  of  r^*, ,  given  by: 

iu>  t  -riti 


J[(^;r)  ■=  Jdtelu>te 


i 


zr 


r  ioj 


r  -ico 


(14) 


We  can  resolve  each  velocity  mode  into  its  two  components 
of  polarization: 


v{t)  =  5.Vv^  ' 


(15) 


where  and  are  Cvo  orthogonal  unit  vectors  in  a  plane 

which  is  normal  to  the  direction  of  k  as  shown  in  Fig.  7. 

If  q  is  an  arbitrary  direction,  and  if  q*k=cos6  defines  the 
angle  between  q  and  £,  then  since  Vg.  is  normal  to  k*,  then  we 
can  pick  one  direction  of  polarization  normal  to  both  k  and  q, 
say  and  the°  write 


=  si*e 


(16) 


and  thus 


<l£- ~  5inl0  <j  v^Cco)!1)  . 


(17) 


We  can  write  <  |  v-j>^  (ca)  |  ^>  as  the  Fourier  transform  of  the 
time  correlation  function  of  Vyv^(t)  aB  follows 


(18) 
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The  equal  tine  correlation  function  Is  given  by  the  equi- 
partition  theorem  of  statistical  mechanics: 


<v> 


(19) 


where  p  is  the  equilibrium  mass  density  of  the  liquid.  The 

time  dependence  of  the  transverse  velocity  modes  is  given 

by  the  Navier-Stolces  equation,  which  states  that  they  decay 

T1  2 

at  a  rate  equal  to  -j^k  ,  where  r)  is  the  shear  viscosity. 
Therefore: 


<Vitt'>Vt')>  =  Te  >  (20) 

and 

<|vri(^}|i>  =  t  (21) 

where  we  have  Inserted  the  wavenumber  and  frequency  dependent 
value  of  the  shear  viscosity  for  r|. 

If  we  Insert  Eqs.  (13),  (17)  and  (21)  into  (12),  we 
obtain 


xr  =  ,7a  I  <i sf.r >[1  - c* )*>i k^).  (», 


In  the  above  equation  we  essentially  have  the  product  of 
two  Lorentzians,  if  TVj*(u>)  has  a  weak  frequency  dependence. 
Since  the  first  Lorentzlan  has  a  much  smaller  width  than  the 
second  one>  we  can  replace  the  second  one  by  its  variation  in 
the  neighborhood  of  o,  that  1b  we  make  the  approximation 
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(23) 


and,  therefore: 


(24) 


B.  Kawasaki  Theory 

Let  us  now  investigate  the  important  special  case,  where 
the  shear  viscosity  is  independent  of  frequency  and  wavenumber, 
and  we  can  write 


V1"’’  1 


(25) 


Then 


ftf^Ppr)  -M^rM  -  P i  .<»> 


and 


<27) 


f 


If  we  assume  an  Ornstein-Zernicke  form: 


<|Sf-«r|l>  -  ,4fr)‘ 


(28) 


then  Eq.  (27)  becomes: 


2T  Ti  fjf  f  dM'-/*1) 
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(29) 


The  second  term  in  the  brackets  of  the  above  expression  does 

not  contribute  because  it  is  an  odd  function  of  y.  To  evaluate 

the  above  integral,  we  substitute: 

.a.  "V 


W  =  I  ~ 

yu.-*  (i-wV/z 

.  aw 
*  "  20^7* 


(30) 


then 


(  1  W  d.W _ 

if  =  1  (.-w)'/4(k^^w),/i 


**i  (i 


Tltx*«'fr  ^  ^  ‘  > 


(31) 


47 


XOhTZ  72-208 


which  is  just  the  Kawasaki  result,  wher«£  K(x)  is  the 

13  14 

Kawasaki  scaling  function  *  ,  given  by: 


(32) 


C.  First  Order  Effect  of  Critical  Viscosity 

Now  let  us  return  to  the  general  case  of  Eq .  (24),  and 
assuae  that  we  can  writ 


V(tj)  =  1- + 


(33) 


where  Ar)^(b))  is  the  critical  wavenumber  and  frequency  dependent 
viscosity,  and  r)Q  la  the  background,  or  ideal  viscosity.  If 
the  anomalous  viscosity  is  a  small  part  of  the  total  vircosity, 
then  to  first  order  we  can  write: 

V"1"  ’  Lv  Ak(<o)]  ’*  %' l'-  ^— ]  •  <«) 


We  will  now  consider  the  wavenumber  dependent  conducti¬ 
vity  at  the  critical  temperature  where  K=o.  Thus  to  lowest 
order 

<*)  iZT  ,  v 


2T 


(35) 


and  the  first  order  correction  is,  then  by  Eq .  (24): 

I  V  J 

1°~~  % 


^  .  (36) 
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By  Bq.  (39)  of  the  previous  chapter,  tre  have 

(T,F' )£(<*>>  rr  +  fy.)  , 


where 


r~  x-  =  r  , 


wr.r) »  V‘*>  •  < 

Therefore,  Eq.  (35)  becomes 


=1t?]  Wdt'e' 

-•*»  -m  -90 


,  io>(t+t) -£ltl  y  It  'I 

e  e  e 


jdtdt  S(t+t')e 


-rjtt 

)*  '  *> 


=  -  7^7  '  <40) 


so  that  Eq.  (39)  becomes: 


a>? = *r.r>  r 


1.‘A 


AlV<rfr>  > 
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where  An*pcr^._^)  Is  Che  Fourier-Laplao.e  transform  of  the 
anomalous  shear  viscosity,  evaluated  at  the  Laplace  frequency 
equal  to  the  relaxation  rate  of  the  q-k  concentration  mode. 
Previously,  we  have  found  that 

2*fi 


(43) 


At  the  critical  temperature  T£  we  can  use  the  formulae: 


<isfr  '>  =  if' 


C-r>’ 


and 


P  =  _E_  if.ri3 

f-e  “>1  1 


(44) 


(45) 


Inserting  Eqs.  (43),  (44)  and  (45)  into  (42)  we  obtain 

2T  n  I  V  t  K  _  ,-n-l—  'll  H.!.\ 

■*4:r"<nr,lV.«,,P'  tU3  )\  > 


where 


r'  =  f 


(47) 


If  we  choose  the  polar  axis  to  be  along  the  direction  of.q, 
and  convert  the  sum  in  Eq.  (46)  to  an  integral,  and  change  to 
dimensionless  variables,  according  to  the  transformation: 


A. 

u=  % 

£t  =  y  =  (uVl-  2UyU)/Z  - 


(48) 
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then 


AXr"~ 


r  ‘rfi'i 


Only  the  first  tern  in  the  brackets  of  the  above  equation 

has  any  dependence  on  q.  When  Integrated,  the  other  terms  will 

give  constants.  Therefore,  let  us  first  examine  the  coefficient 
2kD 

of  the  In  -  tern: 

<1 


*  J  [iT^y.  **""£?*] 


>  .1  '4 

=  ji/iO-/1;  [t  +  ^ vj  ~  *  Trfsn  e»d©  - 

~  •  a  r 


4 1  ' 


Therefore,  we  can  write  the  first  order  correction  to  the 


conductivity  as  follows: 

ZT  $ 


_  _  _1L  JL  i 
f  ”  Hi’  ,51tA  1-  rn T  ^  ^ 


i . 


where  Is  just  a  numerical  constant  given  by 

cx  =  )-  . 

0  _•»  M 

Thus,  to  within  a  constant  Eq.  (51)  is  the  result  we  would 


have  obtained  if  we  had  inserted 


1M-"  A1x 
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into  the  Kawasaki  fornula  for  q»ic: 

h  =  •  <54: 

We  will  evaluate  thiough  the  use  of  two  successive 
transformations.  First  we  will  transform  the  integral  over  y 
to  an  integral  over  u'  by 


udu 


I- A*  -  l  - 


U*+  I +  u'ti*-2u1-2u'z-  2u1u‘l 


i[z(u\u,2)-i  -cu-u’y] 


so  that  Eq.  (52)  becomes 


to  CH*'J  3 

?/v>  ,  C56, 


where  the  region  of  Integration  is  the  diagonal  strip  shown 
in  Fig.  6.  Noting  that  the  term  involving  the  viscosity 
scaling  function  depends  only  on  the  ratio  of  u'  to  u,  we 
transform  to  polar  coordinates  so  that  the  unavoidable  numerical 
integration  be  a  single  integral  Instead  of  a  double  integral. 


We  let 


U  ~  VC.OS  4 

U'~  V  Sin  4 

dudu  -  vdvdfa 
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then 


V/*> 


p  _L  LA£ -  1 2/2- 1  -  24>lfc(.'4tani4> )-  ^Qn{ca^4  *■  4 )*Cri  V ] 

H  irMi.a^CfsVJ  ^  <■  ' 


Jii*^  CpS^J  ' 
O  V*) 


where  is  the  intersection  of  the  ray  <J>  ®  constant 

with  the  line  u’+u  -  1 ,  and  v2(<fr)  is  the  intersection  of  the 
same  ray  with  the  appropriate  one  of  the  two  lines  |u'-u|  *  1, 


V, (.<#)“  (CoS  Sin  41  )  j 


V;Ct)  ~  |ce,s^-  Sin  <fil 


We  can  rewrite  Eq.  (58)  as  follows 
71,1 

cr  ’ (61> 

where  ¥j($)  and  (<*>)  are  the  radial  integrals: 

+  )  sj  %  l*''-  '  -  VW2+]  ,  (62) 


v/dO 


r  l 

^C*)M  ~7 

*  -Wj 


vW2+]  . 


We  will  first  evaluate  4^(40  : 

4W.f  *1^- 7,-v«sz^] 

Vtd») 


(63) 
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L  *.'■*)*] 

"  x  *n  1  cos<f  -S.n«t»  I  H  J 

-7I  s**^?*?  s:  ..»*  f 

2  L  (Cos4>t-S^i4)tCCos4'-Sin<#>)U  J 

=  2£  I52±i^*  1  -  legist* 4 

*  Kn  \  Co*4- \ 

'  i\_6, |ta.  ) I  -  *"  z<#>]  . 


To  evaluate  4»2 (<|»)  ,  we  will  first  establish  a  few  integral 


identities,  namely: 


„  /**v»  r  m*Ox  (n+*)x 

(y*  (  (o^ojc  !  xe  _  £ I 

J  V%,v  <*v  =  J  >e  <<*  =  L'T^y  I"-)1  J. 

*,  £»¥,  ‘r"v* 

•iirv-rM-ar-n ,  ■•» 


which  is  true  fer  n  +  -1,  and 
u. 


J  y  C^  =  JXe4x  =  7  [&»VX-  v>] 


-  4-/w(v,V1)^n(-V-)  . 
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Therefore: 


Vz<*) 


_ —  VCoS22^>  j  CnV  c {v 


/  V*.  \  »  \  Cn  V.  1  +  -L.  - 

^i«(v*v,Kn(-vr )  v/J  v»  -I 

-  i  -  V,Z£«W2]cp522^  +  i  lV*“  V«  JCt,i 

^  Igos^  -*•  Stn<^>  I 

“  —  Cr\  |  Co5^>-Sf»l  i  OJ  I  i 

- 4 j(co3<*> -  sin  <*  )*&  /c*s  c  *  /  -fa&4+*.**i&la>s4  r’*SlW^} 

-  J  (cos4>  ♦  St*  4  )*  —  «»3  <£  -  SwK  <t>  )  \ 

■  l^ltoS^-Sm^/  Atlc^t+su.*! 

2  l  (CoS<f- Sin  4>  i2-  (C«» s4>-t'**<f>)  J 


(tos^-Sin^)1  Ccx>s4-^Sin4> 


rl  (C®S4- -  Sin<i02(Co5<*  +  S.V,  4  J2- 


.  _  /  gc>s<£  ->-  S<>?<^  J 

~  -£*  lctt<t>-Si*<f'l  I  cvs<t  -  si*4  I 
■h  2sr?^>0os^>  £tolCvi<{)  +  SuKt>  (-*-  2Sin<$  CcS<ft&*  jOt>s4'-S<>r>4'  / 

S">2<*j  £n  /Cos  2*/  .  (67) 
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He  note  that  the  Integrand  In  Eq.  (61)  has  a  logarithmic 
singularity  at  ♦  ■  ^  .  Therefore,  we  will  subtract  out  the 
singular  part  of  the  function,  and  integrate  the  remainder 
numerically,  while  the  singular  part  can  be  integrated 
analytically.  In  other  words,  if  we  write 

,v‘ 

CX=J  §(^)d<p  (68) 

a 

where 

£  fen.* CoS3* ) '/ V,(4> ) [<r,M  twV ) - iUeoS •  ,(69) 

and  let  $(<|>)  be  some  function  that  can  be  analytically 
Integrated  such  that 


then  we  write 


Vi  _  r~ 

"  A 


where  the  first  Integration  is  performed  numerically,  and 
the  second  analytically. 

To  construct  i(4>),  we  expand  the  singular  parts  of  $(4>) 

7T 

about  the  singular  point,  4*  *  ■£,  and  evaluate  the  non- 

7T 

singular  terms  at  4>  *  Thus  using  Eq.  (69),  we  can  write 

7T 

Around  <j>  ■  we  can  write 
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Therefore , 


Co S  ♦  Si n  ^ 
cos  -  Sin  ^ 


—  Sin 


2<*J 


(64) 


becomes 


(74) 


and 


(67) 


becomes 


~  1*2  +  (i+£*2)£*l<t-1!zl  +  14-  h  / . 


(75) 


Noting  that 


fWj-WW)*--’®5, 


(76) 


we  can  insert  Eqs.  (76),  (74)  and  (75)  into  (72)  to  obtain 

<|(<*)  -  ^®{2(.<8ls)[l4f»l4-'5l]-l.«3l  ♦  I.LIil&H-lil  ■* 
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*  j.23S-<?  -.7 ifi  -  IxZl<t>-  ?/j  # 


To  integrate  i($)  ,  ve  not'*  that 


o 


=  ({»$- l)f  *-0-2V(.)f  , 


(77) 


(78) 


(79) 


”h  tt/h 

«  4j&?*x^X  =o?[x^x-2x^;X-^xJo 

.  (80) 


Therefore,  by  Eq.  (77) 


% 

/  §C4>M4>  ~ 

o  J 

«  f./-/.3£?/j  =  -0.«7I4.  , 


(81) 


A  graph  of  $($)  -  $(4>)  vs  <j)  is  shown  in  Fig.  8. 
Numerical  Integration  gives  (see  Appendix) 
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-  O.  3  7"V  7  . 


(82) 


Therefore , 


Cx  -  -1. 25-  , 


(83) 


and 


AA.; 


£T 


%  /5n  1i'u"  / 

Thus,  to  first  order,  this  is  equivalent  to  writing 

A?  -.<*<* 


(84) 


r  #<v'* 


(85) 


where 


.cff 


If  -  %  + 


■"  g  " 

Experimentally,  it  has  been  found  that  Eq .  (85)  is  well 
satisfied  if 


(86) 


tf 


(87) 


where 


x 

v 


/.2?  . 


Thus,  if  we  replace  r)^,  by  T)  in  Eq.  (86)  and  solve  for  r| , 
we  find  that 


(88) 


x 

7° 


(89) 
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To  check  the  above  result  with  the  linewidth  measurements 

20 

of  Chang,  Sengers,  Keyes,  and  Alley  ,  It  is  first  necessary 
to  determine  the  value  of  the  parameter  kD ,  the  high  wavenumber 
cutoff.  This  is  done  by  using  the  measurements  of  the  tempera¬ 
ture  dependent  viscosity  of  Stein,  Allegra  and  Allen^.  Both 
investigated  the  binary  mixture  3-methylpentane-nitroethane. 

The  critical  viscosity  was  found  to  have  a  temperature 
dependence  of  the  form 


=  a  U£  +  b  , 


where 


e  - 


7--  Tc 


and  the  best  values  of  a  and  b  were 

d  -  - 1-  Sli  *  IO  <j/r»/cw  -sec. 

b  -  -4-7  «  / o **  ^nn/cm-  sec 

with  a  background  viscosity 


(90) 


(91) 


(92) 


r[0  —  3. £  *  /o^  Q">/c/»-5ec  .  (93) 

The  linewidth  measurements  produced  a  fitted  effective  value 
for  the  viscosity 


1 


-  3  / 

V.  $  *  /O  Jem-  SfrC 


(94) 


Assuming  that  the  correlation  length  obeys  a  power  law  of 
the  form 
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(95) 


then  we  predict  a  temperature  dependence  for  the  critical 
viscosity  as  follows: 

=  5^1-faete)-*333] 


-  £ 


(96) 


where 


b'  =  (97) 

The  linewidth  measurements  were  best  fitted  by  the  values 

•=.  2.56  A 

—  ,(pl(o 

Using  these  values  we  get 


a'^  -  -see.  ,  (" 

or  an  11%  difference  from  the  value  given  by  Stein,  Allegra  and 
Allen. 

If  we  assume  that  our  predicted  value  of  b'  is  the  same  as 
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the  value  by  given  by  the  viscosity  data,  then 


/S’?* 


InLhEj  =  0-333-0.493; 


-  -/-  /5"£ 


*  0./23A  7  . 


aoo) 


The  wavenumber  used  in  the  linevidth  data  was  dependent 
on  the  scattering  angle  of  the  laser  light,  and  was  given  by 


_  *tnn  ,  6 

%~  A  5  l 


(101) 


where  X  *  6328A  was  t^e  wavelength  of  the  laser  light  used, 
and  n  ■  1.38  was  the  index  of  refraction  of  the  liquid.  Thus 


£  =  (jl:M  *  to*  fi~‘ )  Sin  f- 


(102) 


% 


Sin  | 


Using  the  fact  that 


(103) 


JU  U8  --  M.H f 6 


(104) 


and  inserting  it  into  Eq.  (89)  we  obtain 


«-■=  I  1- £,[-<■&*  «*<> ] 


=  {.«28  +  .UH  t  j  . 


(105) 
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The  scattering  angle  ranged  fro*  30°  to  130s  leading  to 


(i) 

Wo  A 


=  1.32. 


max 


(106) 


-  f. 2  2  ,  (io>) 

nt>  n 

whereas,  the  observed  value  is 

/.2?j  <108) 

which,  within  experimental  uncertainty,  is  in  very  good  agreement 
with  theory,  considering  that  it  is  only  a  first  order  theory. 


D.  Viscosity-Dif fusion  Self  Consistency 

In  this  section  we  will  demonstrate  that  self-consistency 
of  the  viscosity  and  the  diffusion  coefficient  demands  that 
no  power  law  divergen:  2  occur  in  the  critical  viscosity.  We 
will  do  this  by  examining  the  Four ier-Laplace  transforms  of 
both  the  viscosity  and  the  diffusion  coefficient.  We  can 
adequately  express  the  Fourier— Laplace  transform  of  the 
diffusion  coefficient  as 


__ 


(109) 


where  we  neglect  the  frequency  dependence  of  the  susceptibility. 

This  is  justified  because  we  know  that  in  the  critical  region 

the  decay  rate  of  the  concentration  fluctuations  becomes 

vanishingly  small,  while  the  decay  rate  of  the  velocity 

fluctuations  remains  finite,  and  possibly  diverges.  To  study 

the  Laplace  transform  of  the  diffusion  coefficient  we  put  a 
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-pt2i 

factor  of  e  into  Sq.  (S)  and  change  the  limits  of 

integration  to  (o ,«)  ,  and  substitute  it  into  Eq.  (109), 
to  obtain 


)><-£ ) £ vtft, )) , 


(110) 


where  we  neglect  the  time  dependence  of  the  concentration 
fluctuations.  If  we  Fourier  transform  in  space 


then 


(111) 


(112) 


Similarly,  by  Eq.  (39)  of  the  previous  chapter  we  have 


(113) 


(114) 


Let  us  now  define  two  characteristic  frequencies 
associated  with  the  critical  Rayleigh  linewidth  and  viscosity, 
pp  and  p^,  respectively.  They  are  obtained  by  replacing  all 
wavenumbers  by  k,  the  reciprocal  of  the  correlation  length, 
as  follows: 


fr  -  KWk' 


(115) 
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(116) 


and  we  know  that  p^»pp. 

We  assume  that  the  critical  behavior  of  the  shear  viscosity 
disappears  at  frequencies  larger  than  p^.  Thus  if  we  look 
at  the  frequency  denominator  in  Eq.  (113),  we  see  that  the 
most  important  contribution  is  from  wavenumbers  smaller  than 
tc ,  so  that  we  can  replace  all  lengths  by  K  \  and  the 
hydrodynamic  value  of  D  is  of  the  order  of  magnitude 


DtflfV)-"  •  <117) 

However,  from  the  frequency  denominator  of  Eq.  (114)  we  see 

that  the  critical  frequency  behavior  of  An  cuts  off  at  Pp , 

a  frequency  much  smaller  than  p  .  Thus  0*(p  )  can  no  longer 

m  k  n 

have  any  critical  behavior.  Therefore,  all  the  power  law 

behavior  on  the  left  side  of  Eq .  (117)  is  taken  up  by  D,  and 

as  a  result  the  shear  viscosity  cannot  have  any  sore  drastic 

behavior  than  a  logarithmic  divergence,  as  was  previously 

22  13 

demonstrated  by  Swift  ,  and  Kawasaki 

Equation  (117)  is  a  general  equation  for  all  fluids.  That 
is,  it  can  be  used  to  determine  the  critical  variation  of  the 
transport  coefficients  J  and  n*  The  liquid-gas  phase  transition 
is  very  similar  to  the  binary  liquid  phase  transition,  with  the 
density  difference  between  the  liquid  and  gas  phases  replacing 
the  concentration  fluctuation  as  the  order  parameter.  The 
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transport  equations  are  also  similar,  with  the  particle 

diffusion  coefficient  D,  replaced  by  che  quantity  —  in 

P  P 


the  liquid-gas,  where  X  is  the  heat  conduct! /ity,  p  the  density. 


and  C  the  specific  heat  at  constant  pressure.  However 
P 

— ,  like  D,  vanishes  proportionately  to  ic.  In  fact,  it  is 
PCP 


possible  to  use  the  sane  arguments  leading  up  to  Eq.  (117)  to 

show  that  the  shear  viscosity  cannot  have  a  power  law  divergence, 

only  a  logarithmic  divergence  at  worst,  and  that  the  decay 

rate  of  the  order  parameter  varies  as  the  third  power  of  the 

21 

wavenumber,  a  result  previously  found  by  Kadanoff  and  Swift 

The  situation  is  more  complex  in  the  superfluid  phase 
18 

transition  .  However,  to  the  extent  of  the  validity  of  mean 
field  theory  and  dynamical  scaling,  it  is  still  possible  to 
use  Eq.  (117)  to  predict  the  behavior  of  the  critical  shear 
viscosity.  In  the  superfluid  we  have  a  complex  order  parameter, 
with  the  superfluid  density  related  to  the  square  of  the  ampli¬ 
tude  of  the  order  parameter,  and  the  super fluid  velocity  propor¬ 
tional  to  the  gradient  of  the  phase  of  the  order  parameter. 

The  diffusion  coefficient  D  of  liquid  helium  in  a  porous 
medium  leads  to  the  anomalous  increase  in  the  thermal  conducti¬ 
vity  in  the  high  temperature  region,  and  describes  the  anoma¬ 
lous  attenuation  of  second  sound  in  the  low  temperature  region. 

-1/2 

It  can  be  shown  the  D  diverges  as  K  ,  so  that  Eq.  (117) 

3/2 

would  predict  that  the  shear  viscosity  vanishes  as  k  in  the 
hydrodynamic  regime,  and  varies  as  the  wavenumber  k  raised  the 
3/2  power  for  k>>K,  if  we  assume  the  validity  of  dynamical 
scaling . 
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In  this  chapter  we  have  shown  how  the  critical  shear 
viscosity  decreases  the  * alue  of  the  relaxation  rate  of  the 
order  parameter  predicted  by  the  aode-coupling  theory  of 
Kawasaki.  He  have  also  explored  the  self  consistent  nature  of 
the  critical  diffusion  and  viscosity*  and  its  extension  to 
other  fluids. 
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IV.  PERTURBATION  THEORY 


In  this  chapter  ve  will  explore  the  effect  of  higher 
order  pertubations  on  the  critical  shear  viscosity  and 
Rayleigh  linevidth.  We  will  present  an  equal  time  justifica¬ 
tion  for  the  decoupling  approximation  in  the  calculation  of 
the  anomalous  shear  viscosity.  Finally  we  will  examine  "vertex 
corrections"  resulting  from  time  dependent  perturbation  theory, 
and  show  that  they  are  very  small  effects  in  the  critical 
Rayleigh  linewidth  and  shear  viscosity. 

A.  Static  Perturbation  Theory 

In  static  perturbation  theory  we  employ  the  fourth  order 
term  in  the  Ginzburg-Landau  equation  for  the  free  energy 
density^.  If  denotes  the  total  free  energy  of  the 

quadratic  terms,  then  the  partition  function  can  be  written 


is} 

w 


Because  F^0^  is  a  diagonal  quadratic  form  in  wavenumber  space, 
it  does  not  induce  any  cross-correlations.  Thus,  to  first 
order  in  b  we  can  calculate  correlation  functions  just  as  before, 
except  that  we  must  introduce,  within  the  expectation  values,  a 
factor 
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s W  »  -  *&££  £>.****  V< 

rrK+**r*>c 


Proceeding  in  this  spirit  we  will  examine  the  effect  of 
the  above  factor  on  the  critical  shear  viscosity,  which  by 
Eq.  (30)  of  Chap.  II  can  be  written  as 


tv  , 

d1f  e 


'  V<4,)  V ct' 5  ^f-r- lti)^  ■ 


If  we  insert  Eq.  (2)  into  the  equal  time  value  ( t21-°'  °f 

the  time  integrand  in  Eq.  (3)  then  we  obtain  a  term  proportional 

to 

which  conveniently  splits  up  into  the  product  of  two  sums,  one 
over  wavenumber  l,  and  one  over  wavenumber  l  .  However,  because 
"q1  is  taken  to  be  along  the  z-axis  both  sums  turn  out  to  be 
antisymmetric  in  the  x-direction  and  as  a  result  are  zero  as 
shown  below: 
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Thus  we  see  that  from  the  staodpolnt  of  static  perturbation 
theory,  the  decoupling  approximation  entering  into  critical 
shear  viscosity  is  justified.  It  is  worth  noting  that  the 
decoupling  approximation  is  valid  to  all  orders  of  static 
perturbation  theory,  because  to  all  orders  the  double  sum  in 
Eq.  (4)  will  be  split  into  a  product  of  two  sums,  leading  to 
the  zero  result  of  Eqs.  (5)  and  (6). 

Unfortunately,  the  above  result  does  not  show  up  in  the 
calculation  of  the  critical  bulk  viscosity.  The  bulk  vis¬ 
cosity  £  is  proportional  to  the  space-time  integral  of  the 
correlation  function  of  the  trace  of  the  stress  tensor^. 

If  we  employ  the  canonical  formalism  of  Chap.  II,  then  it  is 
easy  to  see  that  the  trace  of  the  stress  tensor  has  a  term 
proportional  to  the  square  of  the  order  parameter.  Thus  to 
calculate  the  temperature  dependence  of  the  hydrodynamic 
bulk  viscosity  it  is  necessary  to  cope  with  a  term  of  the  form 


where  f  (  j  l  |  ,  |  JL '  |  )  has  no  angular  dependence,  depending  only  on 
the  magnitudes  of  the  wavenumbers  &  and  it'.  We  see  that  we  no 
longer  have  any  angular  dependence  in  the  integral  to  save  us, 
as  in  che  case  of  the  shear  viscosity.  When  we  substitute 
Eq.  (2)  into  (7),  we  find  that  the  equal  time  value  of  the 
integrand  in  Eq.  (7)  becomes  proportional  to 


Z  21  f  oTijro  <1  Sr  r ‘>‘<1 V  I‘  >l 

Hi  r  r. 
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theory  it  is  necessary  to  consider  relevant  non-linear  terms 
in  the  hydrodynamic  equations,  and  how  they  effect  the  decay 
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of  concentration  and  velocity  modes.  We  will  only  consider 
how  a  concentration  and  velocity  mode  yield  another  concen¬ 
tration  mode,  and  how  two  concentration  modes  yield  a 
velocity  mode.  These  vertices  are  characterized  by 

(Ha) 


The  vertex  (11a)  is  easily  derived  by  considering  the 
convective  term  V* (sv)  as  a  perturbation  to  the  time  rate  of 
change  of  the  concentration  fluctuation.  In  wavenumber  space 
the  hydrodynamic  equation  is 

(ti  *  rr ) y * >  =  (12) 

T'  b 

where  we  treat  the  right  hand  side  as  a  perturbation.  We  can 
think  of  the  zeroth  order  value  8g.°^(t)  as  the  solution  to 
the  homogeneous  equation,  and  s^i‘1^(t)  as  the  first  iteration 
of  Eq .  (12),  that  is  for  t>o: 

¥<*>  =  (13) 


fH> .  (u, 
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To  derive  Che  vertex  (lib) ,  we  must  consider  the 
canonical  form  of  the  stress  tensor.  If  J"(r*»t)  is  the  part 
of  the  stress  tensor  at  space  time  point  (r”,  t)  due  to  con¬ 
centration  fluctuations,  then  we  can  write  the  hydrodynamic 
equation  for  the  transverse  velocity  as  follows: 

,  us) 


where  in  Chapter  II  we  showed  that 


•T  -  s  101  _  C  fir 

J‘i  "  Vas,  V 


and 


with 


n^,}-  is1* 

k:i 


9s  .  . 


:i6> 


(17) 


(18) 


In  Eq.  (17)  we  neglect  fourth  order  terms,  and  also  assume 
that  y*o.  Then 

%  ■  Av^lK-Viv]  , 

and 


(19) 


y—S'.  =  y^'ls  s  ]-ass-  -7^'s  s  . 

fe*  *  ^  L  "  J  *  ‘J-J  i  h=!  * 


(20) 


so  that 


— *  **■*'*’  I  /  2  7  \  — > 

v-tr  =  ^ 


(21) 
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where 


k  —  ozy 


Transforming  to  wavenumber  space,  we  find  that  the  q 
component  of  Is 


(' V  *  ^  ~  K  l]  V-r  V 

0  j* 

-  F.lt(f-rJ5rV 

£  *  ■* 

x-^pr(^  -  y  -  sr  sf-r  •  <23) 

Since  we  are  only  dealing  with  transverse  velocity  components. 
It  is  necessary  to  subtract  out  the  longitudinal  components, 
that  is  the  components  of  (7*T)  along  the  direction  of  the 
wavevector  <f,  so  that 

<*F£-  &UT~i$W*t~t}vow-  <24) 


Using  Eqs.  (15)  and  (24)  we  can  write  the  hydrodynamic  equation 
for  the  transverse  velocity  modes  as 

~  *&Ue~  (25) 


The  solution  to  the  homogeneous  equation  is  just 
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-*  (oi.t  . 

vf  ft) 


.JLj'ir 

— >  fo)  a  js 

(o)  e  '  * 


(26) 


for  t>o,  and  the  first  iteration  is 


r»>* 


s£u-isg.te> 


f* 


(27) 


Having  evaluated  the  vertices  (11a)  and  (lib)  we  are 
in  a  position  to  examine  the  vertex  corrections  to  the 
critical  Rayleigh  linewidth  and  shear  viscosity. 


C.  Vertex  Corrections 

To  examine  the  wavenumber  and  frequency  dependence  of 
the  critical  Rayleigh  linewidth,  we  write  the  fluctuation 
dissipation  formula 

rr<u)’Tk.lM^t<f  :YD)>  -  (28) 

X 

where 


75sZV°f  . 

6  7T 


(29) 


so  that 


(30) 


Thus  we  must  correlate  a  concentration  and  a  velocity  mode 
at  time  t  with  a  concentration  and  velocity  mode  at  time  0, 
as  shown  below: 
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(31) 


Just  connecting  the  lines  to  form  the  bubble 


(32) 


14 

corresponds  to  the  decoupling  mode  approximation  of  Ferrell  , 

13 

discussed  In  Chapter  III,  Section  B,  and  yields  the  Kawasaki 
result . 

To  calculate  the  vertex  correction  we  trace  back  the 
velocity  mode  at  time  t  and  presume  it  was  caused  by  two  con¬ 
centration  modes,  and  trace  back  the  concentration  mode  at 
time  t  and  presume  it  was  caused  by  a  velocity  and  a  concentra¬ 
tion  mode.  We  will  also  assume  that  the  concentration  mode 
was  formed  first,  so  that  the  intermediate  state  consists  of 
three  concentration  modes,  and  no  velocity  modes,  otherwise 
we  obtain  a  correction  of  a  smaller  order  of  magnitude,  as  we 
will  later  show.  Therefore,  we  need  to  calculate  the  following 
diagram: 


t  >-t 


(33) 


We  insert  Eqs.  (14)  and  (27)  into  (30)  and  obtain 


«i 


76 


NOLTR  72-208 


X 


\  * . +~  i  .  i 


(34) 


To  evaluate  the  above  equation  we  take  the  coupling  of  the 
concentration  inodes  shown  by  the  lines  above  the  equation 
and  double  It,  since  by  symmetry  it  must  give  a  contribution 
equal  to  the  coupling  denoted  by  the  lines  underneath  the 
equation.  Therefore,  we  presume  that  as  far  as  Eq .  (34)  is 
concerned 


-fkr-r’t' 


2T  X?xT-r  *>e  e 

*  V;-r 


(35) 


and 


<*  «*»  “  — jpTfF  J 

*  r e  *■  ■ 


(36) 


Note  that 


*T-  -  S'  l*P-",'T1I 


(37) 


The  time  integrals  then  become 
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X  £ 


-  i  itr-rfi  ' 


"...  b* - 


*Uti~)e 


=  |>.f •  (38> 


In  the  above  integrals  -c  neglected  concentration  decay  rates 
with  respect  to  velocity  decay  rates.  Therefore,  if  there 
had  been  any  velocity  modes  propagating  in  the  intermediate 
interval,  then  all  thr je  frequency  denominators  would  contain 
velocity  decay  rates,  and  since  the  velocity  decay  rate  is 
proportional  to  the  square  of  the  wavenumber,  while  the  concen¬ 
tration  decay  rate  is  proportional  to  the  cube  of  the  wavenumber, 
the  sum  in  Eq.  (34)  would  have  an  extra  wavenumber  factor, 
thus  reducing  the  order  of  magnitude.  If  we  substitute  Eqs. 
(35)-(38)  into  (34),  we  obtain 
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£*  £  * 

fa>*e{t*rr] . 


There  ii  another  contribution  to  the  vertex  correction 
It  is  derived  by  tracing  back  the  velocity  mode  at  time  t, 
and  assuming  it  was  formed  by  two  concentration  modes  at  time 
and  further  assuming  that  one  of  the  concentration  modes  was 
formed  by  a  concentration  and  a  velocity  mode  at  time  t" . 
is,  if 


n,«-  <W*’w  -  1 


4 


Vn  f-r-ru'>  - 


then  we  split  up  one  of  the  concentration  modes  at  t': 

If  we  substitute  Eqs.  (40)  and  (41)  into  (30)  we  obtain 
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(46) 


Therefore  the  total  vertex  correction  to  the  decay  rate  of  the 
concentration  nodes  Is  given  by 


e!%>  -  «.  rf’c) 


.  far ) 
a 


-f 


L-i™  far!  I'i0}  f  • 


(47) 


2  3 

The  above  result  is  the  same  as  that  found  by  Kawasaki  . 
According  to  him  the  above  amounts  only  to  a  2.4%  correction 
if  q>>K,  and  0.4%  if  q<<K,  both  in  the  zero  frequency  limit. 
Actually,  it  is  fairly  simple  to  see  that  the  vertex  correction 
to  the  Rayleigh  linewidth  Is  very  small  compared  to  the  zeroth 
order  value 


.-7  (*), 

't  ( 


- 1 


(48) 


because  we  can  write 


Preceding  page  blank 
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fw-7 kz^lt-  ,4F1"KJ  V->  •  '*•> 

where  we  have  multiplied  the  integrand  of  Eq.  (48)  by  the 
correction  function  V^.  jJ’aj) ,  which  by  Eq.  (47)  is  given  by 

V°  Pr^u  V  *r  Jly-rzrTjpry - J 

I^T^Tprj'J l'*’* jV- r J  l"*'1  V* * ^ J • 

r  o  *■  (50) 


V  -  ^(w)  is  a  completely  dimensionless  quantity.  It  is  also 
not  very  difficult  to  see  that  it  is  very  small  for  all 
ranges  of  k  and  q.  If  we  examine  the  integrand  we  see  that 

4 

it  varies  as  k /l  if  2>>k,  in  both  the  very  small  and  very 
large  region  of  q.  If  l  is  very  small,  then  the  integrand 
varies  as  the  inverse  cube  of  the  greater  of  q  and  k.  Xf  l 
and  k  have  the  same  magnitude  or  the  same  direction,  the  inte¬ 
grand  is  zero.  As  a  result  we  expect  (u>)  to  be  of  smaller 

q 

order  of  magnitude  than  (ou )  ,  although  it  has  essentially 

the  same  wavenumber  and  frequency  dependence. 

We  can  use  the  same  techniques  to  examine  the  vertex  cor¬ 
rections  to  the  critical  shear  viscosity.  By  Eq.  (30)  of 
Chap.  II,  the  critical  part  of  the  shear  viscosity  can  be 
written  as 


fc  it 
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Here  we  have  two  concentration  modes  at  time  t  correlated  with 
two  concentration  modes  at  time  0..  as  shown  below 


Joining  the  lines  to  form  the  bubble 


(52) 


(53) 


corresponds  to  the  decoupled  mode  approximation  of  Chapter  II. 
We  will  examine  the  vertex  correction  caused  by  a  velocity 
mode  propagating  In  an  intermediate  state,  as  illustrated  by 
the  diagram: 


(5*) 


To  calculate  the  above  diagram,  we  assume  that  both  concentra¬ 
tion  modes  at  time  t  were  produced  by  a  concentration  and 
velocity  mode  at  some  previous  time.  Therefore,  we  let 


(55) 


and 


r 


a*  o 


(56) 


and  substitute  into  Eq.  (51)  to  obtain  the  vertex  correction 


Atf  ^  “tS ?  IZZ  Jdt'e 
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pITa1 


It  tr  $” 


F(t?-c*Xf'D(f  r)i 

itr-ry  J 


«r)  -l&rfctt-) 

(57) 


Using  the  decoupling  approximation,  and  setting  q  along  the 
z-a  is  as  in  Chapter  II,  we  obtain 


Vt£-r)Cf-eHr-r>  1 
ar-ry  J 


r 


(58) 


Recalling  that  the  uncorrected  critical  shear  viscosity  can 
be  written  as 

&jr  M  ^  Z (fc&t WpXgp [-M f  'r ■*■  ^  r?  >»  (59) 

V 

we  rewrite  Eq.  (58)  in  the  form  of  Eq.  (59)  with  a  vertex 
factor  in  the  summand  as  we  did  with  the  Rayleigh  linewidth, 
as  follows: 


where  W—  <u>)  is  the  dimensionless  quantity  given  by 
q  ,  k 
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Iff  »  -  -for-  o-]7 


(61) 


Because  of  the  proliferation  of  angular  factors  we  expect 
that  the  above  term  should  be  very  small,  although,  it  may 
introduce  aln(kp/q)  term,  since  it  is  necessary  to  cut  off 
the  Integration  over  £  in  order  to  insure  convergence.  As 
in  the  case  of  the  critical  Rayleigh  linewidth,  the  vertex 
correction  to  the  shear  viscosity  has  essentially  the  same  wave- 
number  and  frequency  dependence  as  the  uncorrected  value,  but 
reduced  greatly  enough,  so  that  its  effect  is  negligible. 

In  this  chapter,  we  have  justified  the  use  of  the  de¬ 
coupling  approximation  for  the  shear  viscosity  and  diffusion 
coefficient.  We  have  found  that  static  perturbation  theory 
has  absolutely  no  effect  on  the  critical  shear  viscosity, 
even  though  it  produces  a  oivergent  series  for  the  bulk  viscos¬ 
ity.  After  formulating  the  relevant  techniques  for  time 
dependent  perturbation  theory,  we  examined  the  effect  of 
vertex  corrections  on  the  shear  viscosity  and  Rayleigh  linewidth, 
and  showed  it  to  be  negligible. 
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V.  CONCLUSION 


A.  Summary  and  Discussion 

The  fluctuation-dissipation  theorem  was  a  very  useful 
tool  in  obtaining  the  results  of  this  paper.  It  related  trans¬ 
port  coefficients,  quantities  which  are  proportional  to  the 
dissipation  rate  of  fluctuating  internal  modes  of  the  system, 
to  the  correlation  function  of  the  strengths  of  those  fluctua¬ 
tions.  It  is  intuitively  clear  that  since  a  fluctuation  of 
any  internal  quantity  is  a  deviation  from  the  equilibrium  con¬ 
figuration  of  the  system,  the  stability  criterion,  or  the 
principle  of  minimization  of  free  energy,  demands  a  return  to 
equilibrium.  Therefore,  the  larger  the  fluctuation,  the 

greater  the  restoring  force.  One  manifestation  of  the  fluctua- 

8  9  24 

tion  dissipation  theorem  is  the  Kubo  formula  *  ’  ,  which 

equates  the  transport  coefficient  to  the  time  and  space  integral 
of  the  correlation  function  of  the  relevant  flux.  In  both  the 
shear  viscosity  and  diffusion  coefficient,  the  flux  was  found 
to  be  the  second  order  product  of  two  fluctuating  quantities. 
Therefore,  in  the  calculation  of  the  correlation  functions  it 

14 

was  necessary  to  use  a  decoupling,  or  factoring,  approximation 
which  transformed  a  correlation  function  involving  four  fluctua¬ 
ting  quantities  into  the  product  of  two  correlation  functions, 
each  involving  two  fluctuating  quantities.  The  explicit 
equations  for  tne  correlation  functions  were  derived  from  mean 
field  theory  and  linear  hydrodynamics. 
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14 

These  methods  were  £lrst  used  by  Ferrell  in  calculating 

the  critical  temperature  dependence  of  the  diffusion  coef fx- 

13 

clent.  He  verified  the  result  of  Fawasaki  ,  namely 

& ~  •  (i) 

where 

>"f  ■*  ,  (2) 

is  Stokes'  formula  for  the  mobility  of  a  sphere  of  radius  £ 
moving  through  a  liquid  of  viscosity  r) .  This  result  suggests 
an  intuitive  picture  of  the  phase  transition,  as  consisting 
of  a  binary  liquid  with  regions  of  large  concentrations  of 
solute.  These  regions  can  be  thought  of  as  rigid  spheres 
having  radii  equal  to  the  correlation  length.  As  the  critical 
temperature  is  approached,  the  correlation  length  grows;  that 
is  the  size  of  the  inhomogeneities  grows  larger,  corresponding 
to  a  decrease  in  the  mobility  of  these  regions,  thus  accounting 
for  the  sluggishness  of  the  system,  and  its  lack  of  haste  in 
returning  to  equilibrium.  When  the  correlation  length  grows 
larger  than  the  wavelength  of  the  concentration  fluctuations, 
i.e.  when  q>>£  the  above  formulae  are  no  longer  valid, 
because  D  becomes  a  function  of  q  corresponding  to  "non  local" 
diffusion.  This  change  is  carried  out  by  the  substitution 

>  (3) 
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where  a  Is  sone  numerical  constant  of  the  order  of  unity. 

18  19 

Eq.  (3)  illustrates  the  principle  of  dynamical  scaling  ’  , 

that  all  temperature  dependence  be  expressed  in  terms  of  the 
correlation  length  £»  and  as  T-*-Tc  all  factors  of  £  be  replaced 
by  the  wavelength  of  the  fluctuations.  In  the  case  of  the 
diffusion  coefficient  a  *  3ir/8,  so  that  in  the  limit  T-*-Tc  we 
obtain 

rr 3  *  -  "Sf  •  (4) 

All  of  the  above  equations  assume  a  constant  value  of 
the  shear  viscosity.  However,  recent  experiments  have  clearly 
established  a  critical  temperature  dependence  in  the  hydro- 
dynamic  shear  viscosity  q.  It  was  the  purpose  of  this  paper 
to  theoretically  derive  the  temperature,  wavenumber,  and 
frequency  dependence  of  the  excess  shear  viscosity,  and  deter¬ 
mine  how  it  affects  the  critical  diffusion  in  coefficient. 

In  other  words,  what  value  of  n  should  be  used  in  Eq.  (A)? 

A  Kobo  type  formula  relating  the  shear  viscosity  to  the 
correlation  function  of  the  off-diagonal  component  of  the 
stress  tensor"^  was  employed.  A  canonical  formalism  identify¬ 
ing  the  Ginzburg-Landau  free  energy  density  with  a  Lagrangian 
density  established  that  the  off-diagonal  component  of  the 
stress  tensor  due  to  concentration  fluctuations  is  proportional 
to  the  product  of  two  orthogonal  components  of  the  gradient 
of  the  concentration.  To  evaluate  the  correlation  function. 
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it  was  necessary  to  employ  a  decoupling  approximation,  whose 
validity  was  later  confirmed  by  perturbation  theory.  The 
rest  of  the  calculation  was  rather  straightforward  and  pro¬ 
duced  a  logarithmic  temperature  divergence'1"*.  By  the  princi¬ 
ple  of  dynamical  scaling  a  logarithmic  wavenumber  and  frequency 
dependence  was  predicted  and  actually  calculated,  and  for  the 

shear  viscosity  the  factor  a  in  Eq.  (3)  was  found  to  be  .497, 

25 

confirming  the  general  rule  of  thumb  ,  that  when  factoriza¬ 
tion  introduces  the  correlation  length  twice  (via  the  equal 
time  correlation  function)  ,  a  is  close  to 

The  above  scaling  factor  is  not  directly  valid  in  the 
diffusion  problem  because  it  is  only  applicable  in  the  static 
limit.  Therefore,  it  is  not  correct  to  blindly  insert  the 
wavenumber  dependent  viscosity  into  Eq.  (4)  in  determining 
the  critical  Rayleigh  linewidth.  The  diffusion  process 
depends  upon  the  relaxation  of  current  fluctuations,  which  are 
the  product  of  concentration  and  velocity  fluctuations.  The 
finite  concentration  relaxation  rate  forces  us  to  study  the 
viscous  damping  of  velocity  fluctuations  at  non-zero  frequen¬ 
cies.  This  results  in  a  kind  of  "retardation"  correction  to 


the  effective  critical  viscosity  used  in  the  calculation  of 
the  critical  diffusion  coefficient.  It  is  manifested  by  the 


increase  in  the  scaling  constant  a  ,c  to  be  used  in  the  shear 

eft 

viscosity  which  is  to  be  inserted  into  Eq.  (4).  The  result  is 
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The  self-consistent  nature  of  the  viscosity  and  diffusion 
was  then  studied.  It  was  found  that  the  product  of  the 
diffusion  coefficient  and  excess  shear  viscosity  should  have 
the  sane  temperature  behavior  as  the  inverse  of  the  correlation 
length.  Making  use  of  the  fact  that  the  characteristic  fre¬ 
quency  associated  with  viscous  damping  is  much  larger  than 
the  characteristic  frequency  associated  with  diffusion,  it 

was  deduced  that  the  critical  shear  viscosity  cannot  have  a 

2122 

temperature  divergence  more  drastic  than  a  logarithmic  one 
The  self-consistency  was  then  applied  to  other  fluids,  namely 
the  gas,  and  the  superfluid. 

Finally,  we  applied  the  results  of  perturbation  theory  to 
test  the  reliability  of  the  previous  results  on  the  critical 
shear  viscosity  and  diffusion  coefficients.  The  fourth  order 
term  of  the  Ginzburg-Landau  equation  serves  as  the  perturbation 
in  the  static  case;  it  is  found  to  have  no  effect  on  the  shear 
viscosity,  even  though  it  produces  a  divergent  perturbation 
series  in  the  case  of  the  bulk  viscosity.  Non-lin^ir  hydro¬ 
dynamics  serves  as  the  basis  for  time  dependent  perturbation 
theory,  which  was  used  to  study  the  effect  of  vertex  corrections 
on  the  critical  diffusion  and  shear  viscosity.  Although,  the 
vertex  corrections  had  similar  temperature,  wavenumber  and 
frequency  behavior  as  the  uncorrected  quantities,  their  sizes 
were  diminished  by  at  least  an  order  of  magnitude. 
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B.  Areas  of  Further  Study 

At  first  glance,  it  seems  that  the  techniques  developed 
here  vould  be  suitable  for  studying  the  critical  bulk  viscos¬ 
ity  and  sound  attenuation  In  the  binary  liquid  phase  transi¬ 
tion.  Unfortunately,  It  Is  difficult  to  treat  the  correlation 
function  of  four  concentration  fluctuations  because  of  the 
failure  of  the  decoupling  approximation.  Therefore,  it  will 
be  important  to  learn  how  to  treat  higher  order  correlation 
functions,  and  apply  the  results  to  the  bulk  viscosity. 

A  more  promising  area  of  research  is  the  non-Lorentzian 
behavior  of  the  frequency  spectrum  of  the  concentration  modes. 
This  behavior  will  be  manifested  in  the  frequency  dependence 
of  the  diffusion  coefficient.  The  most  important  contribution 
to  the  frequency  dependence  will  probably  arise  from  the  fre¬ 
quency  dependence  of  the  excess  shear  viscosity,  which  we  have 
already  derived  here. 
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APPENDIX 


In  this  appendix,  the  computer  programs  that  were  used 
for  numerical  calculations  are  presented.  They  were  run  on 
the  INTERCOM  system  of  the  CDC  6400  computer  at  the  U.S.  Naval 
Ordnance  Laboratory. 

Program  FLV283  determined  the  value  of  the  viscosity 
scaling  function  at  zero  frequency,  given  by  Eq.  (109)  in 
Chap.  II. 

Program  SLP283  determined  the  Initial  slope  (minus  1/3) 
of  the  viscosity  scaling  function,  given  by  Eq .  (110)  in 
Chap.  II. 

Program  SCL283  was  used  to  plot  the  viscosity  scaling 
function  for  various  values  of  P,  given  by  Eq.  (100)  in 
Chap.  II.  Simpson's  Rule  was  used  in  the  numerical  integration. 

Program  DIF283  was  used  to  plot  $($)-  5(<f>)  for  various  <f>. 
$(40  is  defined  by  Eq.  (69)  and  $(40  by  Eq.  (77)  of  Chap.  III. 
This  program  used  the  results  of  SCL283,  and  employed  linear 
interpolation  to  calculate  intermediate  values  of  the  scaling 
function.  By  integrating  $(40~$(40  by  Simpson's  Rule,  the 

program  was  able  to  calculate  C-  . 

A 
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-•1 09038 20E4 01 

.330 

-. 10515181E401 

.335 

-• 10137926E401 

•  340 

97709  037E400 

.345 

-•94130832E400 

•  350 

• *9 063533 6E> 00 

.355 

-.8721 4120E400 

•  360 

-.830  59  527E400 

•  365 

-•80S64576E400 

•  370 

-• 7732289 2E4 00 

•  375 

- • 741 28  645E+00 

•  380 

-•709  7  6  49  2E400 

•  305 

- .678  61 542E4Q0 

•  390 

-.6477931 4E400 

•  39  5 

-•61725706 £♦ 00 

•  400 

-•58  69  69  79E4  00 

•  4C5 

-•55689717E400 

•  410 

- »  52 7 008 2 OF* 00 

•  415 

-  •  49  727  48  6E*  00 

•  420 

-. 46767192E400 

•  425 

-.438  17  68  0E4  00 

•  430 

-•  4GS  769  48 £♦ 00 

•  435 

-.37943233E400 

•  440 

- ■ 3  50 1 500 1 £♦ 00 

•  445 

-•3209  09  33E4  00 

.450 

-•29  1699  1  7E*00 

•  455 

-•26251028E40C 

•  460 

-•23333  52  IE* 00 

•  465 

-•2041 68 1 3 E^OO 

•  470 

-• 17500467E4 00 

•  475 

-•  1  458  418  1  £4  00 

•  460 

-•1 1667763E400 

•  465 

-•8751 1 123E-01 

.490 

-• 58335250E-0 1 

.495 

-.28657729E-0! 

WITH  H  «  .005#  C  «  - •  63 2 606 5 4E* 00 

1  3. 39. 20- STOP 
♦  tfCApY. 
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SLP283  14.44.56.  04/12/72 
10  PROGRAM  SLP28  3C  INPUT, OUTPUT) 

15  DIMENSION  GC50) 

20  DATA  PI,  S3/3*  1  41 5926536,  1.7320508076/ 

30  HYPCX)  =  0. 5*AL0GC C 1 .+X)/C 1 .-X) >/X 
40  CIR<X)  =  CPI/2. -ATANC } ./( S3*X) > >/( S3*X) 

50  FHCX)  =  1.  -  X*X 
60  FCCX)  =  1.  ♦  3. *X*X 

70  BR1CX)  =  FHCX)*CHYPCX)-CIRCX))/C 1 6.*X*X) 

80  BR2CX )  =  CFCCX)*CIRCX)-1.)/C24.*X*X) 

90  FCX)  =  7.5*FHCX)*CER1CX)-PP2CX)) 

100  PRINT  100 

1 10  N  =  20 

120  DO  10  I  =  1  ,N 

1  30  Y  =  FL3ATC  I  J/FL3ATCN) 

1 40  GCI)  =  FCY) 

1 50  PRINT  1 10,  Y,  GCI > 

160  10  CONTINUE 

170  100  FORMATC 1 5X,*X*, 14X,*FCX)*> 

180  110  F0RMATC2C 10X, F10. 5) ) 

1  90  S  =  4.+GC  1  ) 

200  N1  =  N/2  -  1 

2  10  DO  20  J  =  1 ,  N  1 

220  S  =  2.*GC2*J)  +  4.*GC2*J+1)  +  S 
2  30  20  CONTINUE 

240  XINT  =  C  S+GCN) )/C  3»*FL3ATCN) ) 

250  PRINT  120,  XINT 

260  120  F3FMATC/5X, * INTEGRAL  OF  FCX)  =  *,F10.5) 
2  70  STOP 
280  END 
** READY. 

FF-TN 

SLP283  14.46.55.  04/12/7? 

X 

.05000 
. 10000 
. 1 5000 
. 20000 
.25000 
. 30000 
. 35000 
. 40000 
. 45000 
. 50000 
.55000 
. 60000 
. 65000 
. 70000 
. 75000 
.80000 
.85000 
.90000 
.95000 
1.00000 


FCX) 

-.00248 
-.009  69 
-.02099 
-.03540 
-.051 78 
-.06895 
-.08577 
-.10125 
-.11 458 
-. 12512 
-. 1 3236 
-• 1 3595 
-. 1 3561 
-  .  1  J 1  1  3 
-. 1 2231 
-. 10895 
-.09075 
-.06729 
-.03774 
0.00000 


INTEGRAL  OF  FCX)  =  -.07908 

1 4.47. 1 3. STOP 
*  +  KEADY • 
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SO.C83  14*41*  36*  0V1C/7C 
10  PM  GUAM  $CLC«3  < INPUT# BUTPUT) 

CO  A TAWHCT)  ■  0* S*ALB6( C 1 .♦¥>/< 1 . -Y> > 

3C  ®UX>  ■  1*  -  X*X 

40  GJt(X>  ■  6H<  X )*( X*ATANH< I #/X ) -  1  .  ) 

50  A(  X  *  P)  ■  SflRT(X*X*X*P) 

60  B(X>  ■  30RTC3.4X) 

70  6C<X*P>  ■  1.  ♦  C  A(X*P>/B<X>>»*C 

80  GJC<X«  P>  ■  GC<X*P>*(  !•-(  A(X*P)/B<X>>«ATAN(B<X>/A<X*P)>> 

90  DtX*P>  ■  A(X*P>**2  ♦  <X*B<X)>**2 

100  BRC<X#P>-Q4<X>*<GJl(XJ*6JCCX#P>>-*B.*D(X*P>/<  1 5«*A( X* P>**C> 
110  F<X*P>  *  -7.54X4X4BRC(X*P>/D<X*P> 

ICO  6<T*P>  *  PC  TANl  T)  *P)/C8  5<  T>4*B 
130  NC  -  50 
1  40  PI  »  3*141 5926536 
150  H  •  •  25/FLi  AT<  N2> 

140  N  •  NC/C 
1 70  N1  ■  P  -  1 
180  PRINT  100 
lfO  08  30  1  »  1*50 

coo  p  »  io***<ci-eoo/io«> 

CIO  Ti  ■  .25  ♦  H 
CCO  61  ■  G(Pl*TlvP> 

C 30  5  *  4.4G1  -2*/< l*4p> 

C40  w  eo  j  »  i#Ni 

250  Tl  -  *25  ♦  H*rL8ATCC*J> 

C60  Gl  »  6€PI*T1#P> 

C  70  TC  ■  *25  ♦  H*n.0ATCC4J^l> 

C60  62  ■  G< PI  *T2* P) 

C90  S  *  2*  *S1  ♦  4. *62  ♦  5 
30C  CO  CCNTINUE 
310  SI  6  ■  PI6H4S/3* 

320  PRINT  1 10*  P*  SI G 
3  30  30  CONTINUE 

340  100  F3RHATC /COX* *P** 28X* *S1 G( P> •/> 

350  1  10  F0RMAK2C  10X*E15*8>> 

3  60  3T8P 
370  END 

**  READY. 

F  S1G<P> 


. 12589254E-01 

•  1  58  48  9  32E-0 1 

•  19952623E-01 
*251 18864E-01 
•3162P777E-01 

•  398 1071 7E-0 1 
. 501 18723E-01 
.63095734E-01 
.79432823E-01 

•  1 OOOOOGOE+OO 

•  1 2589254E+00 
.  158  4893CE  +  00 
.  1 99  52623E+C0 
. C51 1 88  64E+00 
. 31 622777E+0G 

•  398 1071 7E+00 
.501 18723E+00 
•63095734E+00 
•79432823E+00 
• 10000000E+01 
• 12589254E+01 
•158 48932E+01 
. 19952 623E^0l 
.251 18C  64E+01 

•  31 622777E+0 1 

•  398 1071 7E+01 
. 501  18  723E+0 1 
•63095734E^01 
• 79432823E+01 

•  1 OOOOOOOE+02 
. 12589254E^02 
.  1 58  489  32E^02 
. 19952623E^0C 
.251 1 88  64E+02 
•31 622777E+0C 

•  398 1071 7E+02 
.501  18  723E*02 
»63095734E*02 
. 79432823E^02 
. 1000G000E+03 
. 1 2589254E+03 
.158 48932E+03 
• 19952623E403 
.251  18864E*03 
• 31622777E+03 
. 398 1071 7E^03 
.501 18723E+03 
.63095*»34E^C3 
. 7943282CE+03 
. 1 OOOOOOOE^Q4 

1 4.39.41.ST0P 
♦*REAOY. 

LIST 


-•  629 42^  ^8E+00 
-•62861a  1E+00 
-.627593  7E+00 
-• 62631 4.0E+00 
-• 6247 1 582E+00 
-. 62272C 1 7E+00 
• • 62023444E+00 
- • 6 1 7 1 467  6E+00 
-. 61 332375E+00 
-♦ 6086091 6E+00 
-.60282379E+00 
-•5957661 5E+00 
-•  58721  A;75t^  00 
-•57695128T>00 
-. 56474323E+00 
-. 55039 3C5E+00 
-. 53374507E+00 
-•51471 358E+00 
-• 49  330666E+00 
-. 4G964341E+00 
-.  4439  6236E+00 
-•41 061444E+00 
-. 38804418E+00 
-.35875678s>00 
' .32928187E+00 
- • 30 C 1 3327 E* 00 
-•27177 69 6E* 00 
-. 24460 628E+00 
-•2 189282 6E+00 
-• 19495950E+00 
-•  172831 78E  +  00 
-•  1 52  6008 l E^OO 
13426075E+00 
-.1 1 77569  3E*00 
-• 10300045E+00 
-•898  78 408  E- 01 
-.78265445E-01 
-. 68029  655E- 01 
-. 590C9367E-01 
-.51 ! 66899E-01 
4^290852E-01 
-.38298733E-01 
-.  3308  6989 E- 0 1 
- «  28  56 1 4 1 2E-0 1 
-•  24 6 3 732 6E- 01 
-.2123921  IE-01 
-. 18299635E-01 
-. 15759154E-C1 
13565504E-01 
-.1 1 67257CE-0 1 
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D1F283  I7.31.P3.  03/13/7* 

10  DIM  X(IOO)t  Y ( 1 00  J  #  J4CI00) 

20  LFT  N  «  51 

30  F0R  1  •  1  T0  N 

40  READ  /(I),  Y <I) 

50  NEXT  1 

60  LFT  FI  *  3. 1 41 5926536 

70  PRINT 

80  PRINT  "ANO.F(  RAO/PI  >"#  INTEGRAND" 

90  PRINT 

95  LET  L  *  0 

96  LET  R4  n  L0G(PI/4.) 

97  LET  S4<0)  *  -4.*<  .M59-.  71  41*P4-P4»2)/F1  »2 

98  PRINT  L,  S4(0) 

100  F3R  M  «  0.01  T0  0.50  STEP  0.01 

102  LFT  L  «  L  ♦  l 

104  IF  L  ■  25  G0  T0  230 

106  IF  L  «  SO  G0  TO  240 

111  LET  A  «  F*  PI 

120  LET  XO  «  4*(TANCA>)*3 

130  G0SU8  650 

140  LET  T!  *  C8S<A>  ♦  SIN(A) 

150  LFT  T2  *  ABS(C9S(A)-SIN( A) ) 

160  LET  T3  «  A5S< TAN(A»0.25*P! > ) 

165  LFT  T4  «  APS(C0S(?.*A>) 

170  LFT  Jl  «  2.*(L0P<T3)  -  SIN<?.*A5> 

180  LET  J2  *  L0G<T2>»*-L0GCT1>»P*SIN<2.*A>*L0O(T4) 

190  LFT  0  =  S1N( A)*C0S<A) »3 

200  LET  LI  *  C0S( A) *3  ♦  4.*SIN<A)t3 
210  LET  SI  «  J1*(YO-(1/3)*L0GCL1 ))  -  J 2 
220  LFT  S2  «  S|/(D*Plt2) 

22?  LFT  T5  *  L0PC ABS( A -.25*P 1 > > 

225  LET  S3  «  <.♦< .2859-. 71 41 *T5-T54T5> /P I »2 

227  LFT  S4(L)  *  S2  -  S3 

228  G0  T0  250 

230  LFT  S4(L)  «  0.0 

235  00  T0  250 

240  LFT  S4<L>  =  2.*S4CL-!>  -  S4(L-?) 

245  GS  T0  '“SO 

250  PRINT  M,  S4(L> 

260  NEXT  M 

270  LET  C  ■  SACO)  ♦  4.*S4<1> 

280  F0R  J  *  l  T0  <L-2)/2 

290  LET  C  «  2 . * S 4( ?♦ J )  ♦  4.*S4C2*J4|>  ♦  C 

300  NFXT  J 

310  LET  C  *  .01*PI*<OS4<l>)/3.  -  0.8716 
320  PRINT  "CCLAPBOA)  «  "C 

350  ST0P 

650  IF  XO  >*  X(l)  00  TO  690 

660  LET  P  «  (XO-X<| )>/<X(2)-X(| >} 

670  LET  YO  «  < 1 -P) *Y < 1 )  ♦  P+Y ( 2 ) 

680  G0  T0  760 

690  F0R  K  *  2  TP  N 

695  IF  XO  >*  X(K)  GO  TO  740 

710  LFT  P  *  fXO-XCK-l >>/(X(K >-X(K-l > > 

720  LET  YO  «  CI-P)*YCK-I)  ♦  P*Y(K> 

730  GO  T0  760 

740  NEXT  K 

7  SO  LET  YO  «  -.98nR*X0»  < 64{>) 

760  RETURN 

BOO  DATA  0.0,  -.6326 

805  DATA  .01259, -.  6294,. 01  58  5,-.  6286,  .01995# 6276 
810  DATA  . 02512,-. 6263, .03162, -.6247, .03981, -.6227 
815  DATA  .05012, -.6202, .06310, -.61 7|, .0794?, -.6133 
820  DATA  0.10,  - . 60«6 

825  DATA  •  I  259,  -  •  6028  ,  .  :  58  5,  - .  59  58  ,  .  1  99  5,  - .  58  72 
830  DATA  .2512, -.5770-. 31 62, -.5647, .3981, -.5804 
835  DATA  .501 2, -. 5337, . 6310, -.  51 47. .  7943, 4933 
840  DATA  1. 00,  -.4696 

845  OATA  1  .2589, -.4440,  l  .  58  49, -. /I  66,  1  .o«53, -.  388(‘ 

850  DATA  P. 5119, -.3588, 3.1 623,-. 3293, 3.981 l,«.30ri 
855  OATA  5.01  19, -.2718, 6. 3096#-. 2446,7. 9433, -.2189 
860  DATA  10.000,  -.1950 

865  OATA  12.  589, -.  1728,  1  8.849, -.  1  526,  19.953,  -.  1  343 
870  DATA  25. 1 19, -. 1 1 78,31 . 623, 1030,39.8! 1 , -.08988 
875  DATA  $0 . 1 19, -. 07827, 63 . 096, -. n <803, 79 . 4?3, -. 05904 
880  OATA  100.000,  -.05117 

8PS  DATA  125. 89,-. 04429,  ]  sv:.  49, -.f3830,  199. 53, -.03909 

KVO  DATA  251. 19, -.02856, 31  4. ?3,-.r246.4, 394.1  1,.. 02124 
895  DATA  501 . 19, -.01830# X3f  .96, -.01 576, 79/. 33, -.05 1 67 
900  DATA  moo.n,  -.01167 

999  rwr 

**WFAPV . 
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ANGLF<RAD/PI>  INTF GRAND 


0 

-.1 62133 

.01 

-. 1 65618 

.02 

16R  486 

.03 

-.1 70593 

.04 

-.  1  71  789 

.05 

-.  1  7191  7 

.06 

-. 1 70812 

.07 

-. 168294 

.08 

-.164173 

.09 

-.1 58244 

.  1 

-.1 50283 

.11 

-. 1 40045 

.12 

-.127275 

.13 

-.111  6^79 

.  1 4 

-9.30346F-2 

.15 

-7.1 0009F-2 

.16 

-4. 53502E-2 

.17 

-1 .5R554F-2 

.18 

.01 75R 

.  19 

5. 47647F-2 

.2 

9. 4R846F-2 

.21 

.136435 

.22 

. 1 7581 7 

.23 

.204802 

•  2  4 

.200509 

*25 

0 

.26 

-.252907 

.27 

-.331934 

.28 

-.377598 

.29 

-.407413 

.3 

-.42789! 

.31 

-. 4421 1 6 

.32 

-. 45243 

.33 

-.459  595 

.34 

-.464047 

.35 

-.466359 

.36 

-. 466996 

.37 

-.466243 

.38 

-. 46451 4 

.39 

-.461997 

.  4 

-. 4586 

.  41 

-.45453 

.42 

-. 4501 47 

.  43 

-.445387 

.  44 

-.440381 

.  45 

-.43  5251 

.  46 

-.4301 74 

.  47 

-.42518 

.  48 

-.420376 

.  49 

-.41 5882 

.5 

-.41 1388 

CCLAPBDA) 

=  -1.24R04 

**RFADY. 
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